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REPORT  SUMMARY 


Research  is  under  way  '•n  several  areas  of  collision  phenomena  relevant 
to  molecular  laser  operation.  Due  to  the  diverse  nature  of  these  problems,  a 
multi-faceted  approach  has  been  taken.  This  has  led  to  the  study  of  both 
specific  and  model  molecular  systems  in  order  to  understand  the  physical 
phenomena  involved.  We  believe  the  ensuing  results  have  important  implica- 
tions for  the  development  of  high-power,  efficient  molecular  lasers.  In 
addition,  a better  overall  fundamental  understanding  of  molecular  collisions 
has  been  gained  and  this  should  have  application  to  a number  of  other  relaxa- 
tion phenomena.  The  remainder  of  this  report  consists  of  three  manuscripts 
detailing  specific  aspects  of  this  research  performed  in  the  past  year.  A 
brief  summary  of  these  results  is  given  below. 

1 . Kinetic  Analysis  of  the  CO  Laser  System  with  Updated  Rate  Information 

In  collaboration  with  Professor  E.  Fisher,  Wayne  State  University,  we 
have  undertaken  a thorough  examination  of  the  kinetics  in  the  low-temperature 
CO-He  discharge  system.  This  system  has  been  suggested  as  an  efficient 
source  of  high-power  infrared  laser  radiation.  Our  present  research  focuses 
on  an  analysis  of  the  observed  CO  population  distribution  and  its  implica- 
tions for  collisional  rate  behavior.  This  study  has  made  important  use  of  our 
recently  calculated  CO  rate  constants  which  show  significant  scaling  differ- 
ences with  the  older  semiempirical  models.  The  attached  report  presents  a 
simple  and  useful  correlation  of  the  newly  available  rate  information.  This 
correlation  is  presently  being  used  in  CO  laser  analysis  and  the  preliminary 
conclusion  is  that  the  new  rate  scalings  make  a substantial  improvement  in 
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explaining  the  observed  CO  population.  It  is  our  ultimata  goal  to  turn  this 
information  around  and  suggest  design  improvements  for  operating  laser 
systems.  Further  research  in  this  area  will  be  pursued  in  the  coming  year. 

2.  Model  Studies  to  Determine  the  Key  Molecular 
parameters  Controlling  Collision  Behavior 

A variety  of  model  calculations  have  been  performed  to  bring  out  key 
physical  effects  entering  into  molecular  collisions.  A prime  concern  is  the 
determination  of  temperature  and  quantum  number  scaling  behavior.  This  infor- 
mation is  very  important  in  laser  analysis  and  other  general  relaxation 
phenomena.  A new  and  more  accurate  computer  code  has  been  written  to  replace 
our  earlier  version  fj.  Chem.  Phys.  64,  5291  (1976)].,  Calculations  are 
presently  under  way  and  will  be  report  d at  a later  date. 

3 • Stochastic  Theory  for  Molecular  Collisions 

In  the  past  few  years  the  development  of  powerful  effective  Hamiltonian 
methods  has  greatly  simplified  molecular  collision  calculations.  Indeed  many 
previously  impossible  problems  are  now  manageable  with  these  methods  and  we 
have  applied  them  in  the  research  of  this  contract.  However,  the  present 
available  theory  has  limitations,  particularly  concerning  the  size  of  problem 
that  can  be  handled.  Th, s restriction  is  very  important  in  many  practical 
applications  such  as  lasers.  Therefore,  we  have  undertaken  a new  approach 
based  on  stochastic  theory  to  handle  large  (many-quantum  level  ) collision 
problems.  The  attached  manuscript  presents  the  stochastic  theory  formuLation 
for  vibration-rotation  collision  state  changes.  Tne  paper  discusses  the 
physicaL  content  oi'  the  new  simplified  theory  and  suggests  numerical  methods 
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for  its  implementation.  Applications  are  under  way  and  will  be  pursued 
further  in  the  coming  year. 

1 ' 

! 

! Vibration-Rotation  Collisions  in  H, -H2 

The  H2-H2  system  provides  a valuable  prototype  example  for  many 
collisional  effects  of  importance  in  lasing  systems.  Detailed  three- 
dimensional  quantum  mechanical  calculations  were  performed  on  this  system 
with  emphasis  on  the  exchange  of  rotational  and  vibrational  quanta  between 
the  molecules.  Collisional  flux  maps  were  studied  to  determine  the  referred 
pathways  for  vibration-rotation  transitions.  This  study  has  led  to  valuable 
insight  into  the  collision  of  identical  molecules.  The  details  of  this  work 
are  included  in  the  attached  manuscript.  It  vrould  be  valuable  to  extend  this 
study  to  the  collision  of  non-identical  molecules  such  as  D2-H2,  since  such 
situations  commonly  arise  in  many  applications  such  as  lasers.  Further  con- 
sideration of  this  work  will  he  treated  in  the  coming  year. 
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ABSTRACT 

A new  correlation  of  the  recent  theoretical  CO-He  vibration-transla- 
tion rates  by  Vert er  and  Rabitz  is  presented.  Using  a functional  form 
suggested  by  the  SSH' theory,  a simple  analytical  correlation  function  is 
obtained  for  3000°K  > T > 100°K  and  quantum  state  change  Ai  = 1,2,3.  In 
addition,  a semi -empirical  formula  for  Ai  = 1 rates  is  obtained  by  fine 
tuning  the  correlating  formula  to  agree  with  the  available  experimental 
1 to  0 rates.  In  this  fashion  optimal  use  was  made  of  the  best  available 
theoretical  ar.d  experimental  information.  The  result  is  compared  with 
rates  estimated  by  other  procedures. 
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CORRELATION  OF  VIBRATION-TRANSLATION 


RATES  FOR  CO-He 
by 

S.  H.  Lam*  and  H.  Rabitz** 


I.  INTRODUCTION 

In  the  study  of  vibrational  kinetics  of  CO  laser  systems,  quantitative 

information  on  various  vibrational  rates  are  required.  A reexamination  of 

the  CO  V-T  (vibration-translation)  information  at  this  time  is  appropriate 

because  of  its  importance  in  laser  studies  and  because  new  theoretical^  and 

2 

experimental  results  have  recently  become  available.  In  the  present  paper, 
we  shall  be  solely  concerned  with  V-T  rates  for  CO-He  collisions,  particu- 
larly at  low  temperatures  T and  high  vibrational  quantum  numbers  i.  Our 
objective  here  is  to  construct  and  present  a simple  semi -empirical  formula 
for  the  estimation  of  these  rates  based  on  the  best  theoretical  and  experi- 
mental information  at  our  disposal  at  this  time.  The  proposed  semi-empirical 
formula,  in  addition  to  being  of  practical  usefulness,  exhibits  certain 
general  features  where  are  interesting  theoretically  and  are  expected  to  be 
present  in  other  similar  systems. 
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We  shall  denote  k^(i,T;Ai)  to  be  the  de-excitation  rate  of  a CO 
i'.r.lecule  from  i to  i-Ai  at  temperature  T.  The  corresponding  Ai  to  0 
;ates  shall  be  denoted  simply  by  k^T  g(T). 

Experimental  data  for  K^g(T;exp}  is  relatively  abundant,  covering 
a wide  range  of  temperature.  ’’  Experimental  data  for  i/1  rates  is 
generally  scarce,  and  for  CO-He  the  totality  of  data  appears  to  consist 

• C 

of  only  5 points  obtained  by  G.  Hancock  and  Smith''  (for  i=9 ,1 0, 11 ,12,13) . 
Thus,  estimation  of  rates  for  i>l  and  Ai>l  in  practical  applications  must 
rely  primarily  on  guidance  from  theoretical  considerations. 


II.  THE  SSH  THEORY  AND  THE  BR  ESTIMATION  PROCEDURE 

The  vibrational  levels  E-j  of  CO  may  be  written  approximately  as 
follows 


Ei  = ike  1 


- 

1 


-xe(i-l) 


i * 0,1,2,., 


where  0y  = 3123°K  is  the  characteristic  vibrational  temperature  and 
-3 

xe  = 5.98  x 10  is  the  anharmonicity  which,  inspite  of  its  smallness, 

exerts  a strong  influence  on  the  V-T  rates. 

The  basic  theory  on  vibrational  rates  was  developed  in  a sequence  of 

6 7 

two  papers  by  Schwartz,  et  ai.  ’ commonly  referred  to  as  the  SSH  theory. 
The  SSH  result  for  Ai=l  V-T  rates  can  be  written  in  the  following  form: 


kVT(i,T;l)  = (TH|  (i)A^T(y,T)  (2.1) 

where  u^(T)  carries  the  dimension  of  k^1  and  contains  mainly  the  total 
collision  rate.  The  function  J^(i)  is  proportional  to  the  square  of  the 
appropriate  oscillator  matrix  element  which,  for  anharmonic  oscillators, 
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can  be  approximated  by 

VU  ’ T7.V  (2.2) 

VI 

and  A-j  (y,T)  is  the  so-called  adiabaticity  factor.  The  parameter  y is 
defined  by 

y * y(i.T)  * 0-2xei)  (2.3) 

where  8 depends  on  molecular  parameters  and  the  assumed  interaction 
potential.  For  a repulsive  exponential  interaction  potential  with 
characteristic  length  i,  0 is  given  by 

0 = 167rV2k0j/h2  ‘ (2.4) 

where  y is  the  reduced  mass  of  the  collision.  Using  the  analytical  result 

o 

of  Jackson  and  Mott  and  approximately  accounting  for  3-dimensionality 

effects7,  the  SSH' theory  yields  the  following  expression  for  A^(y,T): 

AEi 

A^y.T)  = e 2kT  f(y,z)  (2.5) 


where  AE-j  = (l-2x  • )k0v  is  the  energy  defect  and  f(y,z)  is  the  thermal 
averaging  integral 
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In  the  limit  as  z -*■«>,  1(5, y,00)  reduces  to  $inh~2(y/^5) , and  fly,**) 
becomes 


s*co 

f(y>M)  = j e"^sinh‘2(y/^C))dF 

(2.9) 

The  function  f(y,°°)  has  been  numerically  evaluated  by 
who  suggest  the  following  curve-fitting  formula 

g 

Keck  and  Carrier 

f(y>“)  = FKC(y),  0 - y < 20> 

(2. )0) 

FKC(y)  = 1/2  (3-e"2y/3)e"2y/3 

(2.11) 

Eq.  (2.11)  is  commonly  preferred  over  the  steepest-descent  result, 

8(n  y7^)1^”3^3,  (y>>l),  because  of  its  uniform  validity  in  the 
indicated  range. 

An  estimation  procedure  for  V-T  rates,  apparently  first  used  by  Bray10 

for  high  temperature  and  later  more  extensively  developed  by  Rich11  and 
12 

co-workers  , for  low  temperature  proceeds  as  follows.  First,  the  function 
FKC(y)  is  adopted  to  represent  f(y,z).  The  value  of  0 (or  £)  is  either 
estimated  or  taken  from  tabulated  data  (e.g.,  Herzfeld  and  Litovitz  ). 

The  factor  wj(T)  is  then  determined  empirically  by  fitting  tne  resulting 
formula  for  i=l  to  the  relatively  abundant  experimental  data  on  k^Cnexp). 
If  additional  i f 1 data  are  available,  they  are  used  as  a consistency 
check.  We  shall  call  this  method  the  BR  estimation  procedure. 

Note  that  the  BR  estimation  procedure  relies  primarily  on  F^y)  for 
the  major  i and  T scaling,  and  the  final  formula  contains  a single  molec- 
ular parameter,  0.  A large  number  of  kinetic  calculations  have  adopted 
this  procedure  because  of  its  simplicity  and  because  of  the  absence  of 
any  better  estimates. 
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Recently,  Verter  and  Rabitz^  (VR)  performed  extensive  numerical 
calculations  of  V-T  rates  for  CO-He  collisions  over  the  range  3000°K 

1 T > 100°K  and  1 <.  i <44,  including  1 = 1,  2,3.  These  ab^  initio 
calculations  provide  for  the  f<rst  time  detailed  theoretical  information 
hitherto  unavailable.  Although  further  refinements  in  the  calculations 
are  possible  (at  considerable  additional  expense),  these  results  are 
believed  to  be  the  best  theoretical  CO-He  rates  presently  available. 

Verter  and  Rabitz  correlated  their  computed  results  in  the  follow- 
ing form: 

4 

kVT(i ,T;Ai)  = i exp  £ Bn(T,Ai)in_1  (2.12) 

n=G 

To  cover  the  range  3000°  > T > 100°K  in  10Q°K  increments  for  Ai  = 1,  2,  3, 
a total  of  4E0  coefficients  (Bp ' s ) were  needed.  In  this  form  the  VR 
results  are  somewhat  cumbersome  to  use  and  difficult  to  generalize  or 
interpolate.  Comparison  of  VR  and  experimental  k^g(T)  rates  shows  agree- 
ment with  a factor  of  about  5 at  high  temperatures  and  a factor  of  about 

2 at  low  temperatures.  Comparison  between  the  VP  rates  and  the  commonly 
used  BR  rates  shows  substantial  differences,  particularly  at  high  values 
of  i and  low  temperatures.  These  differences  can  have  significant  effects 
on  kinetic  modelling  in  CO  laser  studies  as  shall  be  pointed  out  later. 

III.  NEW  CORRELATION  OF  THE  VR  DATA 

Motivated  by  Eq.  (2.1),  we  propose  the  following  expresion  as  an 
alternative  to  Eq.  (2.12): 

AiAE-j 

kVT(i ,T;Ai)  = u>Ai(i,mAi(i)e  *kT  ^.(yj) 


(3.1) 
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where  H^(i)  remains  as  given  by  Eq.  (2.2)  and  are  chosen  to  be,b 

i2(i)  = fi-j  (i  )^-j  (i -1 ) , 43(i)  = Jl1(i)S,1(i-l)S.1(i-2)  (3.2) 

The  function  F^(y,T)  is  to  be  chosen  under  the  guiding  principle  that 
w^i (i ,T)  should  have  the  minimum  i dependence.  We  shall  retain  y as 
defined  by  Eq.  (2.3)  as  the  primary  correlation  variable.  The  parameter 
6 affects  y only  by  a constant  factor  and  need  not  be  specified  in  the 
correlation  procedure.  In  fact,  we  shall  presently  determine  an  effec- 
tive 6 as  a by-product  of  our  correlations. 

Inspection  of  the  VR  data  appropriately  reduced  and  p*otted  immedi- 
ately indicates  that  the  hi  = 1 high  temperature  data  can  indeed  be 
roughly  correlated  by  FKC(y)  provided  we  choose  0 = 2.64  * 106.  Fixing 

O 

0 at  value  (which  corresponds  tc  l = 0.3  A),  the  following  form 


for  FAi(y,T)  is  proposed. 


FAi<y,T)  ' f 


,-ay*-ctb]  (y-y*) 


1+e  ala^"a2  ) 


where  a,  a-j,  a2,  b-| , b2,  b^,  and  y*  are  curve-fitting  parameters.  The 
first  term  is  suggested  by  FK(.(y)  and  dominates  in  the  y < y*  region. 
The  second  term  dr  •'nates  in  the  y > y*  region.  The  parameters  are 
found  to  be 


a(Ai=l ) = 2/3, 
a1  = (2/3)\ 
y*  = 21.5 
h - fm\Z 

b2  ~ VTV  ’ 


a(Ai=2)  = 1.03,  a(Ai=3)  = 

3 2 *"  5.5 

b1  = C.  525  (1  + ” j 
k - 150 
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Note  that  6 = 2.64  x 106  °K  is  obtained  as  a consequence  of  setting 
a(Ai=l ) = 2/3. 

Using  this  specific  choice  of  F^. (y,T),  the  values  of  u>A.(i ,T)  are 
then  computed  and  typical  distributions  are  shown  in  Fig.  1.  It  is  seen 
that  the  resulting  (o^(i,T)'s  generally  do  exhibit  a weak  i dependence. 

As  shall  be  seen  more  clearly  later  in  Fig.  3,  the  exceptional  cases 
are  <^(1 ,100°K) , u>2(i<10,  T<400°K)V  and  cc-3(i<10,  T<1000°K)  where  the 
correlation  becomes  locally  poor.  It  is  believed  that  for  Ai  = 3 the 
form  for  <L(i)  assumed  in  Eq.  (3,2)  is  mainly  responsible  for  these 
large  local  aberations.  For  the  other  cases,  the  local  aberations  may 
be  reduced  by  fine  tuning  the  parameters  b2  and  b^,  Since  the  original 
VR  data  for  these  cases  is  likely  to  have  the  largest  numerical  errors 
(noise),  such  refinements  were  not  considered  meaningful  or  worthwhile 
at  this  time. 

Fig.  2 shows  Tu^(T)  vs  T where  ^(T)  is  the  averaged  value  of  w ,. 
(10<i<40,T).  Fig.  3 shows  F^(y,T)  as  given  by  Eq.  (3.3)  vs  y as  solid 
lines.  The  original  VR  data  represented  by  F^. (y,T)u^. (i  ,T)/w^(T)  may 
also  be  displayed  on  this  diagram.  However,  because  of  the  scale  of 
the  diagram,  only  the  data  with  (i ,T)/w^^ (T)  substantially  different 
from  unity  will  deviate  from  the  solid  lines  and  is  shown  by  the  dotted  lines. 

IV.  INCORPORATION  OF  EXPERIMENTAL  DATA  FOR  Ai  = 1 

The  weak  "scatter"  of  w-j(i,T)  about  i^(T)  is  comparable  to  the 
scatter  of  the  original  VR  calculations.  Thus  we  may  justifiably  replace 
w-j(i,T)  in  Eq.  (3.1)  by  w, (T)  and  consider  the  resulting  formula  as  an 
alternative  to  Eq.  (2.12)  when  Ai  - 1. 


Borrowing  the  idea  from  the  BR  estimation  procedure,  we  propose  the 


VT 

following  semi-empirical  formula  for  k (i.T,l): 

AEi 

kVT(i  ,T;1  ;emp)  = {T;exp)2.^  (i)e2!sT 


F-j  (y»T) 


(4.1) 


The  function  w-j(T;exp)  is  to  be  determined  by  fitting  the  available 

experimental  data  on  k^Q(T;exp)  while  F^(y,T)  determined  from  the  VR 

VT  -1/3 

theory  is  retained.  Fig.  4 is  a conventional  plot  of  leg  k^  g vs  T , 

3 

showing  the  high  temperature  experimental  data  of  Millikan  and  White 

4 

and  the  moderate  and  low  temperature  data  of  Miller  and  Millikan  and 

2 

the  recent  data  of  Drozdoski , Young,  Bates,  and  Hancock.  Also  shown 
are  the  original  VR  rates. 

Using  the  experimental  data  shown  in  Fig.  4,  the  values  of  w-|(T;exp) 
are  computed  and  are  shown  in  Fig.  5.  Also  shown  are  c^(l,T)  and  w-j(T) 
for  comparison  purposes.  The  u)^(T;exp)  data  can  adequately  be  fitted  by 
• 2 

u-j (T,exp)  = A^q^  , 100°K  < T < 800°K  (4.2a) 

= a(~-V/3  800°K  < T < 3000°K  (4.2b) 

where 

A = 6 x 1013  cm3/sec-molecule 
= 1 x 10"1 °em3/sec-mole. 


In  the  indicated  temperature  range,  Eqs.  (4.1)  and  (4.2)  is  the  proposed 
semi -empirical  CO-He  V-T  rate  formula  wincn  contains  the  most  up-to-date 
theoretical  and  experimental  information  available  at  this  time. 

Fig.  6 shews  kVT(i ,T;1 ;emp)  and  the  corresponding  BR  rates  (using 


0=2.64  x 106K)  for  T = 100°K,  200°K,  and  300°K.  Also  shown  are  the 


five  data  points  of  G.  Hancock  and  Smith^  which  are  believed  to  have 
been  obtained  at  room  temperature.  Unfortunately,  these  data  points  do 
not  discriminate  between  the  two  rate  estimates  which  diverge  from  each 
other  substantially  only  at  higher  vibrational  levels.  Thus,  definitive 
confirmation  of  the  VR  rates  must  await  additional  discriminating  experi- 
mental data  to  become  available. 

v.  discussion 

We  have  obtained  a new  correlation  of  the  theoretical  v'R  rates  for 
CO-He  in  a form  suggested  by  the  SSH  theory.  For  the  Ai  = 1 case,  we 
have  further  proposed  to  replace  u)-|(i,T)  by  u^(T;emp)  given  by  Eq.  (4.2) 
so  that  the  resulting  semi -empirical  formula,  Eq.  (4.1),  reproduces  the 

Vi  , 

experimental  data  on  k-j  0(T;exp). 

Using  harmonic  oscillator  wave  functions,  the  SSH  theory  showed 
that  the  square  of  the  matrix  elements  can  be  written  as  i,j(i)  --  i, 

^(i)  = i(i-l).  Our  choice  of  S^.(i)  here  rests  entirely  on  analogy, 
accounting  for  the  effects  of  anharmonicity  only  through  the  denominator 
of  Eq.  (2.2).  We  believe  that  the  poor  correlation  of  the  low  i data 
for  the  Ai  = 3 case  is  primarily  a reflection  of  the  inadequacy  of  our 
choice  of  &3(i ). 

In  addition  to  the  constant  0 used  in  the  definition  of  y,  the 

function  F^(y,T)  contains  seven  curve-fitting  parameters:  a,  a-j,  a2» 

* 

y , b-j , b2»  b^.  The  values  of  a and  0 are  not  independent.  Ry  setting 
a(l)  = 2/3  and  requiring  F-|(y,T)  to  agree  with  FKr(y)  for  10  < y 20, 
we  deduce  0 = 2.64  x 106  °K,  yielding  from  Eq.  (2.4)  an  effective  inter- 

o 

action  length  % - 0.3  A which  is  found  to  be  in  good  agreement  with  the 


actual  interaction  potential  (in  the  strongly  repulsive  wall  region) 
used  in  the  original  VR  theory.  The  other  curve-fitting  parameters 
are  probably  also  related  to  the  interaction  potential  used,  but  such 
relationships  have  not  been  identified. 


Using  0 = 2.64  * 106  °K  fee  CO-He  collisions,  the  values  of  y for 
the  temperatures  range  considered  covers  the  range  57  > y > 5.  Since 
FKC(y), as  given  by  Keck  and  Carrier,9  is  a valid  approximation  to  f(y,z) 
only  for  y < 20,  we  have  numerically  computed  f(y,z)  as  given  by  Eqs. 
(2.6)  and  (2.7)  Our  results  show  that  for  y < 20, f(y,z)  is  indeed 
accurately  represented  by  FK(,(y).  For  y > 20,  f(y,z)  exhibits  a weak 
T dependence  similar  to  F^(y,T),  but  diverges  rapidly  from  it  with  in- 
creasing y. 

An  interesting  feature  of  the  VR  rates  is  the  distinctly  more  rapid 
variation  of  FA^(y,T)  with  y for  the  first  few  levels  whpn  y > y*  = 21.5. 
Ignoring  the  Ai  = 3 case  which  was  discussed  earlier,  this  distinctive 
"boundary  layer"  structure  for  low  values  of  i can  be  clearly  seen  in 
Fig.  3.  This  behavior  is  accounted  for  by  the  denominator  of  the  second 
cerm  in  Eq.  (3.3),  and  the  parameters  b2 (T ) and  b,(T)  i ere  chosen  to  fit 
the  Ai  = 1,  1 00° K and  200°K  data.  Since  no  such  behavior  is  anticipated 
by  the  SSH  theory,  we  have  carefully  examined  the  possibility  of  syste- 
matic numerical  errors  in  the  original  VR  calculations  and  concluded 
that  it  is  uni  ike .y.  However,  because  of  the  complexities  of  the  numer- 
ical procedures  u'-ed,  this  possibility,  although  remote,  cannot  be 
entirely  ruled  out. 


Except  for  this  boundary  layer  structure  (for  y > y*),  the  VR  data 
is  qualitatively  consistent  with  the  interpretation  that  the  effective 


CVf^iT' 
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value  of  Jt  is  a variable  and  decreasing  function  of  the  "energy  defect" 
of  the  collision. 

While  the  reasonable  agreement  between  the  original  VR  rates  with 
k^0(T;exp)  as  shown  in  Fig.  4 is  encouraging,  it  must  be  emphasized 
that  ro  definitive  verifications  at  high  values  of  i are  yet  available. 

5 

The  consistency  check  with  the  data  of  G.  Hancock  and  Smith  is  incon- 
clusive. Nevertheless,  it  is  clear  that  the  semi-empirical  rates 
(Fqs.  (4.1)  and  (4.2)  proposeo  here  rest  on  firmer  grounds  tnan 
the  8R  rates.  Recently,  Lordi  and  Rich^  measured  CO  population 
distributions  in  a CO-A  laser  system  at  300°K  up  to  i = 35  and  compared 
their  data  with  kinetic  calculations.  It  was  reported  that  calculations 
using  BR  rates  yielded  generally  poor  results,  but  calculations  using 
the  ac[  hO£  assumption  of  replacing  FKC(y)  by  a constant  yielded  much 
more  reasonable  results.  Since  the  VR  rates  exhibit  a generally  weaker 
i dependence  than  .the  BR  rates  (see  Fig.  6),  their  observation  lends 
some  support  to  the  VR  rates.  The  full  implications  of  these  new  rates 
on  CO  laser  calculations  will  be  discussed  in  a later  publication. 
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FIGURE  CAPTIONS 


Fig.  1 


Fig.  2 
Fig.  3 


Fig.  4 
Fig.  5 


Fig.  6 


Typical  distributions  of  u>. . ( i ,T)  as  a function  of  i. 

til 

A flat  distribution  indicates  good  correlation. 


w^(T)  ?s  a function  of  T where  w^.(T)  is  m^OJ) 
averages  over  the  levels. 


The  solid  lines  are  F^(y,T).  The  dotted  lines  are 
FAi(y»T)«Ai(i,T)/wA.(T)  which  represent  the  original 
VR  data.  Because  of  the  scale  of  the  diagram,  only 
the  data  with  w^^(i,T)/u)^^(T)  substantially  different 
from  unity  is  shown. 


Comparison  of  theoretical  and  experimental  rites. 

a)-|(T,exp),  w-|(T)  and  w^(l,T)  as  a function  of  T.  The 
solid  dots  are  the  recent  data  of  Drozdoski,  Young, 
Bates,  and  Hancock.  The  open  circles  are  uj^(T).  The 
crosses  are  t^(l,T).  The  solid  curve  is  the  data  of 
Millikan  and  White.  Straight  line  segments  are  given 
by  Eqs.  (4.2). 


Comparison  of  present  semi-empirical  rates  with  BR 
rates  at  low  temperature.  Note  that  the  BR  rates 
show  a much  stronger  i dependence. 
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ABSTRACT 


Three-dimensional  quantum  mechanical  calculations  in  the  effective 
potential  approximation  have  been  made  on  the  para-hydrogen  system.  At  low 
temperatures  vib -rotational ly  inelastic  collisions  were  examined  while 
breathing  sphere  calculations  were  used  to  probe  the  high-temperature  regime. 
It  was  found  that  simultaneous  vibrational  and  rotational  processes  contribute 
to  the  overall  mechanism  of  vibrational  relaxation.  Both  intra-  and  inter- 
molecular  energy  transfer  is  possible  in  the  present  calculation,  and  the 
corresponding  cross  sections  are  examined  in  detail.  Rates  are  calculated 
from  the  cross  sections  and  compared  with  experiment. 
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I.  INTRODUCTION 

Processes  involving  energy  transfer  in  bimolecular  collisions  are 

frequently  categorized  as  rotation-translation  (R-T),  vibration-translation 

(V-T),  vibration-vibration  (V-V)  and  vibration-rotation-translation  (V-R-T). 

A study  of  a system  of  gaseous  molecules  must,  if  it  is  t.o  be  complete,  focus 

simultaneously  on  all  these  modes  of  energ}  transfer.  It  is  also  most 

desirable  to  treat  the  problem  within  a quantum  mechanical  framework. 

Vibrational  excitation  and  de-excitation  in  collisions  is  an  important  process 

in  gas  laser  systems  and  relaxation  phenomena.*  Research  on  the  contributing 

physical  factors  in  such  energy  transfer  is  desirable.  Understanding  the  role 

of  rotational  inelasticity  accompanying  the  vibrational  inelasticity  is  also 

believed  crucial  in  interpreting  some  relaxation  phenomena.  The  analysis  of 

the  scattering  properties  of  a system  demands  a reasonably  accurate  knowledge 

of  the  interaction  potential  in  order  to  properly  reproduce  the  dynamics.  As 

a practical  matter,  the  computational  expense  must  remain  within  a reasonable 

limit  if  detailed  analysis  is  to  be  possible. 

In  this  context,  the  para-hydrogen  system  is  an  attractive  one  with  a 

view  to  examining  the  processes  mentioned  above.  A study  of  hydrogen  is 

useful  in  itself  and  may  provide  a model  for  more  complex  systems.  The  purely 

vibrational  problem  has  been  approached  in  the  distorted  wave  formalism  by 
2a 

Calvert,  who  calculated  the  ratio  of  de-excitation  rates  for  H2-H2  and  H2-D2. 

A collinear  study  of  near  resonant  V-V  processes  '2-D2  collisions  was  made 
2b 

by  Alexander  using  a variety  of  intermolecular  potentials,  including  the 
one  used  in  this  study.  More  recently  rotations  and  vibrations  have  been 

2c 

included  in  a three-dimensional  semiclass  leal  study  by  Fisher  and  Billing. 

We  have  chosen  to  solve  the  quantum  scattering  equations  for  this  system  in 
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the  effective  potential  (EP)  approximation.^3  With  this  approach  the 

dimensionality  of  the  coupled  equations  that  must  be  solved  is  greatly  reduced 

3b 

in  comparison  to  the  exact-close-coupled  equations.  In  this  fashion  the 

large  numbers  of  channels  involved  in  a vibrotational  calculation  can  be 

4 

manageably  handled.  The  shock  tube  experiments  or  L.'efer  and  Lutz  and  Dove 
and  Teitelbaum"’  as  well  as  the  low- temperature  stimulated  Raman  scattering 
experiments  of  Ducuing  et  al.^  for  H2  allow  a comparison  with  some  of  the 
calculations . 

We  have  examined  the  para-H2 — para-H2  collision  from  two  points  of  view. 
Firstly,  the  molecules  are  treated  as  two  breathing  spheres  (i.e,,  they  are 
restricted  to  the  rotational  state  J = 0).  In  the  second  case  tne  coupled 
vibration-rotation  equations  are  solved  in  the  EP  approximation;  the  effect 
of  including  rotational  sub-levels  on  the  ground  and  first  excited  vibrational 
levels  was  then  studied.  This  latter  work  serves  to  complement  the  partial 
pictures  of  purely  rotational  and  breathing  sphere  calculations. 

Section  II  reviews  the  theoretical  framework  pertinent  to  this  work. 

In  Section  III  the  intermolecular  potential  is  discussed.  Sections  IV  and  V 
examine  the  cross  section  and  rate  behavior  of  the  V-T  and  V-V  and  V-R-T 
processes  as  well  as  the  associated  relaxation  phenomena.  The  final  Section 
VI  draws  conclusions  from  the  study. 
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II.  THEORETICAL  FRAMEWORK 

3 

The  EP  approximation  has  been  extensively  discussed  earlier,  both  with 

7 8 7 9 

regard  to  theoretical  development  and  practical  validity.  The  formulation  ’ 

is  extended  to  the  case  of  two  identical  vib-rotors  in  the  present  treatment. 

The  Hamiltonian  of  this  system  is 


H - T + H0(int)  + V 


(1) 


where  T is  the  kinetic  energy  operator,  V the  interaction  potential,  and 
Ho(int)  the  unperturbed  vib -rotational  Hamiltonian  operator  of  the  colliding 
molecules.  Effective  eigenstates  of  Ho(int)  are  written  as  }n1j,njj2), 
with  n^  describing  the  vibrational  and  rotational  quantum  states  of 
molecule  i*l,2.  The  interaction  potential  V is  usually  expanded  in  a sum 
of  spherical  harmonics  over  the  orientation  of  angles  of  the  two  molecules 

^ ^ O 

V(fi  ,r2,R)  <=  ^ 2.  AA1A2A(rir2R)^imi,e2in2U1t2Am) 

l\&il  raiiU2ra 


(2) 


where  A ^ ^ i are  f-he  expansion  coefficients,  rj  and  r2  define  the 
vibrational  coordinates  of  the  two  molecules  in  the  space-f'.  red  coordinate 

-4 

system,  and  R is  the  vector  joining  the  centers  of  mass  of  the  two  molecules. 
This  expression  is  considerably  simplified  in  the  body-fixed  axes^  shown  in 
Fig.  1.  Further  discussion  of  the  potential  is  contained  in  Section  III  below. 
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The  effective  potential 


10 


matrix  elements  are  given  by 


<nij1n2J2jveff(R)|n1‘ji'nJ’j2'>  « N0  exp(iH) 

I <nin*l  Ajti  ni'n2'>  <i  i JallTje  ^jJI  J i' Ja'>  O) 

txl2j C 1 * 


where 

N0  - ([j,][j2nji']tj2’])'1/4 


■n  «{bi+ j*-  ji’  - u\  + ji*+  j8'+  Ji  + j*}  — 


<JxJaliT 


r[JiHJ*HJi'Hj2,]W]'|l/*  /jl  jl'\ 

(4tt)3  J \0  0 o / 

\o  0 0 / 

[k]  » 2k+l,  and  (II*)  is  a 3-J  symbol.** 

The  expansion  coefficients  have  been  integrated  over  the  vibrational  coordinates 

with  respect  to  the  vibrational  states,  (rjn.)  * <p  (r) 

x n. 


‘‘n»lAi,M<R)lDi",«'>  'll  dr‘dr*  q’n1(r')  'pnJ(r»)  Ai,i2i(r'r«R>  V(Tl)  V(r,) 


(4) 


Equation  (4)  assumes  separability  of  the  molecular  vibration-rotation 
wave functions;  this  approximation  can  easily  be  relaxed  if  necessary. 


In  the  case  of  homonuclear  diatomics  the  additional  symmetry  under 
nuclear  Interchange  requires  that  the  total  wavef unction  be  symmetric.  The 


© ' 
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observed  cross  section,  a0^8,  is  a statistically  weighted  sum  of  symmetrized 
cross  sections 


o = a+W  + a W 


where  the  + and  - signs  refer  to  the  total  symmetry  of  the  spatial  part 
of  the  wavefunction.  For  hydrogen,  which  is  a boson,  these  weights  are  given 


1+1 
21+1  ’ 


2 1+1* 


where  I is  the  total  nuclear  spin  (0  for  para-H2). 

The  wavefunction  is  then  expanded  in  a symmetrized  basis  set. 


k,i,(l)  "W2’ 1 Vcos9>  (6) 


where  0 is  the  polar  angle  of  R.  This  results  in  the  following  form  for 
the  symmetrized  wavefunction 


*±(R’l'2)  ■ T I Vj.-.j.W 


where  the  sum  is  over  all  the  indices.  In  this  basis  set  the  symmetrized 
Ep  matrix  elements  become 
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<Xn1j1n2j2|v+(eff)(R)  liVjjVj,') 


6X<e'  { (1  + 6n!n26  j,  j2)  (1+6n1,n2‘6  jj'jj'  )} 


1/2 


X {<n1j1n2j2lveff(R)|n1'j1'n2‘j2'>  ± (-1)^ Jxna j2| Veff|n2' j2V J,')} 


(8) 


The  coupled  differential  equations  for  the  system  are  then 


. Jd  / _di  _ UtLll  } 


['£{ 


2 ^ !>  dR2  R2  } + En1j1  + En2j2  *"]  UXnt  j.  r.2  j2  ^ 


- £ Unijln2j2|v±efflXnl,j1,n2’jl,)U*n  ,jiVj  ,(R) 

ni'ji'at’J,' 


(9) 


where  p.  is  the  reduced  mass  of  the  colliding  pair. 


E , are  the  internal 

n.  1 . 


energies  of  the  molecules  i=l,2,  and  E is  the  total  energy  at  which  the 

molecules  collide.  The  ground  state  of  each  molecule  Eqo  is  the  zero  point 

energy, 0.268  eV,  and  the  total  energy  satisfies  the  relation  E > 2Eoo.  The 

asymptotic  solution  to  the  Eqs.  (8)  yields  the  scattering  amplitude  from  which 
i 

the  symmetrized  S~  matrix  may  be  obtained 


u"  , (R)  ~ j~6  ,6  ,6,6,  exp  (i  k 

\ 2^  / L n1n1  n2n2  jiji  j 2 J 2*  ' ni 


J »n2  J2 


R cos  0^ 


n,n2'  n2ni’  j 2 J i'  j 1 J 2*  ^ ^ ni 


±6  ,6  , 6 . , , 6 

n n n.r»*  1 4* 


j 1 n2  j 2 


R cos 


(n-  - 9))J 


p exp  (^i  K^,  j rJ  (nt  j tn2  j2  Ej'  j )'n2'  j2'  |q)  ± f(nJj1n2j2->  n2'j  2,nj*  j j'  | tt  - 0 )j 


(10) 
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where  k is  a wave  vector.  The  symmetrized  scattering  amplitude 

jin2 $2 

may  be  defined  as 


£~ (oi  nx' jx’nj' jj’Je ) = f (ni  j:n2j2^  n2' JiV 


± f(n1j1n2j2->n2*j2,n1’j1,{  tt -8) 


( ^^,'/^Cn1  jxn2  j2  n^  j I'n^  j2*  ^ 


» 

x l (21+l)[(u  (-1)1  ‘j.jJCl*  (-D* 

jC-o 


The  total  cross  section  is  then  given  by 


° (nijinaJj^ni'ji'na'ja' ) “ £ 


l 

TT [MI3£I] 

nljln2.J2  _ IjlHjjl  J 


00 

1 («  + ‘)[i*  (-1)*  6ni»a  • ‘..V  63,'iJ 


.V. 

x i6  „,6.  - s’  (n1j1n2jI;nl,j1'n2'j2')| 

1 n^i'  jj*'  n2n2'  j2j2  z 


^'£}g>&£&  fcAyS'J 
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The  above  results  can  readily  be  applied  to  the  case  of  two  breathing  spheres 
by  setting  ji  * j2  * Ji' m it  * 0 and  retaining  only  the  term  ix  - 12  = A = 0 
in  the  potential  in  Eq.  (3). 

After  generating  total  cross  sections  a°^s  as  a function  of  kinetic 

energy,  rates  of  excitation  and  de-excitation  can  be  obtained  as  a function 

12 

of  temperature  by  averaging  over  the  kinetic  energy.  The  appropriate 
excitation  rate  constant  is  given  by  the  expression 


3/2  1/2  * 

WT)  - («-)  i^Y  ‘ e*r(ih-Zi),n)  J (t+Ej-Bj)*  (.+E  )exp(-e/ltT)ae 

0 

(13) 


where  E ^ and  e * E-E^  is  the  final  kinetic  energy.  The  abbreviated 

indices  i,j  label  the  molecular  levels  and  K is  Boltzmann's  constant.  In 

1 

a similar  fashion,  for  de-excitation  (E^  > E^),  the  expression  becomes 


(14) 


Equations  (13)  and  (14)  are  related  by  detailed  balance  and  in  this  paper  the 
* 

units  of  k(T)  are  cm3  sec"1  molecule" 1 .The  various  equations  developed  in 
this  section  will  be  applied  in  Sections  IV,  V and  VI. 


III.  THE  INTERMODULAR  POTENTIAL 


The  interaction  potential  between  two  hydrogen  molecules  has  been 

13 

discussed  in  some  detail  in  the  literature,  and  much  work  is  still 
14 

proceeding.  Unfortunately  for  our  purposes  most  potentials 
developed  thus  far  are  unsuitable  for  application  in  a study  of  vibrations 
since  they  examine  the  H2-H2  surface  keeping  the  intramolecular  H-H  bond 
distance  at  the  equilibrium  value  of  ~ 1.4  Bohr.  Since  the  effective  potential 
matrix  elements  in  Eq.  (3)  must  be  obtained  by  integrating  over  the  vibra- 
tional coordinates  in  Eq.  (4),  the  potential  surface  must  be  known  as  a 

function  of  these  coordinates.  Currently  the  only  available  potentials  that 

13f 

incorporate  this  feature  are  the  Silver -Karp lu3  London  and  Valence-Bond 
'L3st  13 

(VB)  surfaces.  s As  was  demonstrated  elsewhere,  the  inherent  simplicity  of  the 

la'tter  surfaces  precludes  a high  degree  of  accuracy.  However,  these  surfaces 

13  f 

are  believed  to  be  qualitatively  correct,  and  they  should  allow  for  a 
determination  of  the  basic  physical  effects  dominating  the  collision.  This 
goal  ie  in  accord  with  the  purpose  of  this  research,  which  is  the  determination 
of  the  collisional  mechanism  that  is  operative  in  hydrogen  relaxation. 

The  potentials  were  generated  as  a function  of  the  angles  0lf  62,  and 
(cpj  -<p2),  and  the  distances  rlf  r2  and  R.  They  were  then  fit  by  a least 
squares  procedure^  to  the  expansion  in  spherical  harmonics  in  Eq.  (2), 
including  all  terms  up  to  A2i*..  This  expansion  incorporates  the  principal 
components  of  the  potential:  the  spherical  term  A0oo  anc*  the  various  short- 

and  long-range  anisotropies  A20j,  A022  and  A224.  In  computing  the  coupling 
matrix  elements,  the  H2  molecule  was  treated  as  a Morse  oscillator, ^ and  the 
integration  in  Eq.  (4)  performed  numerically. 
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A few  additional  comments  are  necessary  concerning  the  potential.  A p 
the  H2  molecule  is  stretched  to  a bond  distance  greater  than  3 Bohr,  the 
potential  becomes  decidedly  more  anisotropic  than  can  be  adequately  expressed 
by  a least  squares  fit  to  the  six-term  expansion  employed  above.  The  number 
of  vibrational  states  n that  can  be  included  is  therefore  limited  to  those 
whose  wavefunctions  cpR(r)  die  out  sufficiently  by  r = 3 Bohrs.  This  is 
not  a severe  restriction  since  the  bas„.c  trends  are  discernible  within  this 


framework . 
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IV.  COMPUTATION  OF  CROSS  SECTIONS 
A.  General  Comments 

18 

The  Gordon  algorithm  was  used  in  solving  the  coupled  scattering 
equa'.ions,  Eq.  (10).  The  utilization  of  symmetrized  basis  states  in  Eq.  (7) 
requires  solving  fewer  coupled  equations.  This  procedure  also  results  in 
the  calculation  of  a~ 


+ 

j” 


(nl  j ln2i2  '*  nl 


C (nl  j Xn2  J 2 J 3 2*  ) ~ O'C11!  j l n2  j 2 n2'  J 2*n  1*  Jl*  ) (^5) 


which  is  an  inseparable  sum/difference  of  direct  and  exchange  cross  sections 
respectively  (neglecting  the  generally  weak  interference  terras).  Note  that 
for  para-H2  the  statistical  weight,  W = 0,  and  only  matrices  need  be 
calculated  to  obtain  a on  Eq.  (5).  In  actual  computations  the  basis  set 
used  In  calculating  even  or  odd  Jl  valued  matrices  is  different  because 

of  the  nature  of  the  symmetrized  potential  matrix  in  Eq.  (8).  The 
matrices  corresponding  to  even  or  odd  l values  behave  similarly  as  a function 
of  Z,  but  they  can  have  differing  contributions  to  the  total  cross  section. 
Within  r.L'- her  set,  however,  it  is  often  sufficient  to  compute  the  matrices  at 
a uniformly  spaced  grid  (rather  than  each  Z value)  since  the  behavior  is 
smooth  as  a function  of  A. 

The  total  energies  of  the  calculation  were  on  the  range  1.05  eV  < E < 3.12 
for  the  case  of  two  breathing  spheres  and  0.5  eV  < E < 1.35  eV  for  two  vibrotors 
Since  vibrationally  excited  states  are  closed  below  1,05  eV,  only  purely 
rotational  information  is  obtained  in  the  lower  regime.  A more  thorough 
treatment  of  this  aspect  is  to  be  presented  elsewhere,  and  only  brief 
comments  concerning  pure  rotations  will  be  made  below.  Vibrational  transit! ms 
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inherently  lead  to  extremely  small  cross  sections  (save  the  near-resonant 

V-V  cases)  particularly  near  thresholds.  The  difficulties  of  calculation  in 

19  1 

these  cases  have  been  discussed  elsewhere,  and  one  must  use  tight  tolerances. 

Nevertheless,  a certain  degree  of  ’’noise"  still  enters,  and  we  have  found  it 

expedient,  with  the  rrte  calculations  of  the  following  section  in  mind,  to 

20 

smooth  them  using  a spline  fit. 

B.  Breathing  Sphere  Transitions 

The  maximum  basis  set  employed  in  these  calculations,  in  the  notation 
nxn2,  is  showr.  in  Fig, 2 (a)  Tba  actual  basis  set  in  any  particular  calculation 

Hh 

depended  on  the  number  of  states  open  at  that  energy.  The  matrices  were 

tested  for  basis  set  convergence  and  at  varying  tolerances  to  ensure  numerical 
stability.  It  was  generally  necessary  to  include  two  closed  channels  for 
convergence . 

Due  to  differences  in  tie  two  potentials  used,  especially  in  their 

slopes  at  small  intermolecular  separations,  cross  sections  for  the  various 

2b 

inelastic  transitions  can  sometimes  differ  by  an  order  of  magnitude.  The 
steeper  slope  of  the,  VB  potential  gives  rise  to  larger  cross  sections. 

However,  the  behavioral  aspects  of  the  cross  sections  from  the  two  potentials 
are  quite  similar.  In  particular,  transitions  involving  the  change,  of  the 
same  number  of  quanta  are  grouped  together  and  are  generally  separated  from 
other  groups  by  approximately  an  order  of  magnitude.  This  is  illustrated  in 
Fig.3(a)which  shows  the  results  of  calculations  using  the  VB  potential.  Each 
band  encompasses  the  de-excitation  cross  sections  involving  the  same  number 
of  quanta,  An,  exchanged  with  translation 


I 
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H2(n:,0)  + H2(n2,o)  ■*  H2(nx\o)  + H2 (n2’,0) 

An  * (nx  *r  n2 ) - (ni’  + nj') 

The  energy  deficit,  AE  * E _ + E „ - E ,_-E  is  different  for  the 

J nxO  n20  ni'O  n2  0 

various  transitions  within  each  band,  and  this  together  with  the  different 

coupling  effects  accounts  for  the  spread  of  the  band*.  Figure 3(b)  and  all 

following  figures  depict  the  results  of  calculations  using  the  Vb  potential 

since  we  feel  that  this  potential  is  more  in  agreement  qualitatively  and 

quantitatively  with  the  true  physical  potential  (cf.  Section  VI  below). 

It  is  quite  evident  in  Fig.3(a)that  the  near  resonant  V-V  transitions, 
for  which  An  = 0,  have  greatly  enhanced  cross  sections  and  are  well  separated 
from  those  for  which  An  > 0.  The  cross  sections  for  An  > 0 tend  to  overlap 
considerably  in  the  threshold  region.  The  individual  cross  sections  generally 
rise  with  energy  at  different  rates  as  shown  in  Fig.3(b''for  a few  typical 
cases.  This  behavior  can  give  rise  to  different  temperature  scaling  which 
will  be  discussed  later. 

The  detailed  behavior  of  the  cross  section  ^n  Fig .3(b) is  rather  complex, 
but  some  comments  can  be  made.  For  instance,  it  is  apparent  that 
a03~>  11  < °03-»  02  even  thou8h  energy  transferred  is  slightly 

smaller  than  Aj?o3-*02  , Both  these  transitions  result  in  a net  change  of 

An=l.  However,  physically  these  constitute  different  processes.  Neglecting 
the  weak  exchange  term,  the  case  03 ->02  corresponds  to  a change  of  one 
quantum  on  the  second  molecule,  while  03">11  requires  one  molecule  to  gain 
a quantum  and  the  other  to  lose  two  quanta.  This  latter  case  is  clearly  less 
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favorable  due  to  diminished  coupling  effects.  A similar  situation  occurs  in 
the  case  of  00  > cil-»  00  * Caution  is  called  for  in  any  generalization 

of  these  results  since  it  is  apparent  that  ° 22-^11  > °22->02  ’ at  *east  at 
lower  energies. 

C.  Vibration-Rotation  Transitions 

In  the  calculation  of  vibration-rotation  matrices  by  close  coupling 

or  even  effective  potential  methods,  the  principal  deterrent  is  the  number  of 
channels  that  must  be  Included.  While  in  the  case  of  hydrogen  this  problem  is 
reduced  due  to  the  large  rotatioial  spacing  and  weak  coupling,  even  in  this 
case  the  problem  rapidly  becomes  serious.  This  is  obvious  in  Fig.2(b)which 
shows  the  high  density  of  the  rotational  states  involved.  In  tackling  a 
problem  of  this  size,  a certain  judicious  choice  must  be  made  concerning  which 
states  to  include  in  order  to  achieve  convergence  of  the  desired  S*  matrix 

elements.  Another  practical  restriction  is  imposed  through  the  time  required 

18  21 
per  calculational  step,  which  is  a function  of  the  number  of  states.  In 

practice  any  basis  larger  than  20-25  states  imposes  excessive  demands  in  terms 

of  computer  time.  Calculation  of  each  matrix  in  the  basis  of  Fig. 2(b) 

* 

took  8-10  minutes  in  double  precision  on  an  IBM  360/95  computer  since  the 
18 

tolerances  had  to  be  consistently  maintained  at  5 x 10"6. 

Tests  of  convergence  were  performed  by  varying  the  states  in  the  basis. 
Purely  rotational  transition  matrix  elements  generally  converged  on 

inclusion  o*  all  states  up  to  the  energetically  highest  open  level.  Coupling 
decreases  rapidly  for  large  multi-quanta  transitions,  giving  rise  to  extremely 
small  cross  sections  for  transitions  like  0606-»0000.  Inclusion  of  such  very 
high  0j!0j2  states  did  not  affect  the  principal  rotational  and  vibrational 


trans it ions. 
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The  actual  basis  set,  therefore,  consisted  of  the  twenty  states  shown  in 
Fig. 2(b) 

Figure  4 (a) shows  several  purely  rotational  cross  sections  ae  a function 
of  energy.  The  principal  ones,  for  which  Aj  = 2,  lie  mainly  in  the  region 

, >2  o 2 

10  2 A < a < 10  A , and  Aj  > 2 cross  sections  decrease  with  increasing  Aj. 

Figure 4{b)is  a profile  of  de-excitation  rotor  transitions  at  a total  energy 

of  1.06  eV.  It  is  apparent  that  for  inelastic  rotational  cross  sections 

(c  . the  largest  member  corresponds  to  (ji  + j2)  - Oi'  + j2')  = 2« 

' nlJln2J2  nl  Jl  n2  J2  ' 

In  addition,  since  the  lower  quantum  states  are  energetically  closer  together, 
one  readily  sees  for  ji  < j2  that 

> aj!j2  -»  ji  j2  " 2 

Mn.ny  of  those  trends  were  previously  observed  in  atom-molecule  collisions 

such  ns  He-F2  and  1°  rigid-rotor  studies  on  H2-H2.  Although  the  trends 

9a 

are  the  same,  the  magnitudes  of  the  previous  H2-H2  cross  sections  based  on 
a different  potential  are  somewhat  smaller.  The  same  general  behavior  is 
seen  for  taose  transitions  that  involve  one  vibrationally  excited  molecule, 
with  an  additional  interesting  aspect  that  can  be  seen  in  the  comparison  of 
the  two  collisions 


tl2  (02 ) + H2(00)  -*  li2(00)  + H2  (02 ) 


H2  (12)  + H2(00)  -*  H2(10)  + H2  (02) 


(a) 

(b) 


*'*/.*£■*  ' ' ’ ' - . 
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Case  (a)  is  really  an  elastic  collision  due  to  the  enormous  contribution  of 
the  exchange  term  o (0200 -» 0200 ) in  Eq.  (15).  In  case  (b)  the  exchange  term 
is  extremely  small  since  the  process  would  be  a vibrational  exchange. 

Therefore,  case  (b)  oust  rely  on  the  direct  term  cr  (1200 1002)  which  itself 
describes  an  exchange  of  rotational  quanta  between  the  two  molecules.  The 
net  result  is  that  case  (a)  has  a much  larger  cross  section  than  case  (b). 

Figure  5 represents  a typical  profile  of  vibrotor  de-excitation  cross 
sections  at  a total  energy  of  1.2  eV.  The  transitions  to  the  left  of  state 
0100  on  the  abscissa  are  vibrotor  cross  sections,  while  pure  rotor  cross 
sections  (neglecting  the  very  weak  exchange  process)  are  to  the  right.  It  is 
directly  visible  that  as  the  rotational  quantum  change  Aj  increases  beyond  2, 
the  cross  sections  decrease  rapidly.  This  emphasizes  the  importance  of 
coupling  effects  and  shows  the  small  magnitude  of  near  resonant  vibrotational 
sections  which  are  necessarily  accompanied  by  large  changes  of  rotational 
quanta . 

As  can  be  expected,  the  largest  cross  sections  are  those  for  which 

1m!  is  0 or  2.  The  largest  vibration-relation  (although  not  substantially) 

cross  sections  from  0010  and  0212  are  to  0000  and  0202  respectively.  Similarly 

the.  largest  cross  sections  from  0210  and  0012  are  to  0202  and  0000,  respectively. 

These  latter  cases  involve  simultaneous  vibrational  and  rotational  transitions, 

while  the  former,  purely  vibrational.  Some  parallel  observations  on  the  He-H2 

21 

system  have  been  made  recently. 

It  must,  however,  be  noted  that  the  total  energy  of  1.2  eV  in  Fig.  5 
is  the  same  for  all  the  cross  sections.  Hence  the  various  cross  sections  are  at 
different  kinetic  energies  and  this  accounts  for  the  generally  low  magnitude 
of  transitions  from  level  0414  which  is  barely  open  at  1.2  eV. 
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For  the  comparable  breathing  sphere  transition  a is  slightly 

01  ^ 00 

VR 

smaller  than  a v as  has  been  observed  in  the  He-H*  system  by  Zarur 

0100-^0000  ■■  J 1 

and  Rabitz.  Recent  work  using  the  atom-breath? ng  sphere  model  for  tht 
22  23 

He-CO  and  other  model  systems  has  shown  the  approximation  to  give 
reasonable  results  under  certain  circumstances.  Additional  comments  on  this 


matter  will  be  made  in  Section  VI. 
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V.  KATE  AND  RELAXATION  BEHAVIOR 
A.  Breathing  Sphere  Rates 

The  rates  were  calculated  from  the  cross  sections  by  Eqs.  (12)  and  (13) 
over  a large  range  of  temperatures.  The  rates  are  presented  in  Fig.  6 
analogous  to  the  cross  sections  in  Fig. 3 (a).  The  behavior  qualitatively 
follows  that  which  can  be  deduced  from  the  cross  sections.  However,  it  must 
be  noted  that  differing  threshold  behavior  gives  rise  to  slightly  different 
temperature  dependence.  At  high  temperatures  the  dependence  is  Landau-Teller 
while  non-linear  behavior  on  a log  k(T)  vs.  T”1^3  plot  shows  up  below 
- 600°  K. 

The  greatly  enhanced  rates  for  the  near  resonant  V-V  processes  are  a 
consequence  of  their  extremely  large  cross  sections  in  comparison  to  the 
An  > 0 transitions.  The  V3  potential  indicates  a difference  be  ween  the 
rates  for  V-V  and  V-T  processes  of  approximately  a factor  of  103-104.  This 
is  a very  significant  difference  and  it  is  not  clear  if  the  inclusion  of 
rotational  states  would  tend  to  diminish  it.^c’^  Nevertheless,  it  is  evident 
that  ladder  climbing  mechanisms  (i.e. , transitions  of  the  type 
njn?  ->  nj  -l,n2  + l)  are  extremely  important  in  the  rapid  and  efficient 
transfer  of  population  to  higher  vibrational  levels. 

Experimental  measurements  of  the  rate  of  self  relaxation  of  vibrationally 
excited  H2  molecules  by  the  simulated  Raman  scattering  (SRS)  method  find  no 
difference  in  the  rates  for  ortho-or  para-H2  above  300° K Below  room 

temperature  the  para-H2  - para-H2  rate  is  slightly  smaller  than  the  ortho  case. 
Figure  7 presents  a comparison  of  the  results  of  this  calculation  with  both  the 
SRS  experiment  (para-H2)  on  the  range  50°-500° K and  the  shock  tube  data^'“* 
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from  1000° -3000° K (normal  H2).  Of  the  two  potentials  used  the  results  of 

the  VB  potential  are  closer  to  the  experimental  ones, which  suggests 

that  this  is  the  more  accurate  one.  The  VB  potential  has  been  used  in  a 

2b 

collinear  study  of  D2-D2  near  resonant  V-V  processes  by  Alexander.  In 
comparison  wic'tt  several  ab  Initio  potentials,  it  was  observed  that  the  VB 
potential  gave  smaller  cross  sections  and  this  was  attributed  to  the  slow 
decrease  of  the  potential  at  long  range.  While  the  agreement  with  experiment  at 
high  temperatures  is  substantially  better  than  at  low  temperature,  some 
discrepancy  still  exists.  However,  the  breathing  sphere  rate  behavior  is  at 
least  qualitatively  accurate. 


B.  Vibration-Rotation  Relaxation 

Rates  were  generated  in  a similar  fashion  for  V-R-T  transitions  below 
300°  K.  The  restriction  to  this  range  was  necessitated  by  the  availability 
of  only  low  energy  cross  sections.  The  results  of  these  calculations  are 
plotted  in  Fig.  7 in  the  following  manner.  First  a rotationally  summed 
de-excitation  rate  is  defined  by  the  relation 


tot 

klJiOj2 


(T) 


£ kUiOj2^Oj1'Oj2*(T) 

U'h' 


(16) 


These  rates  are  shown  in  the  inset  of  Fig.  7.  The  fact  that  rates  from  the 

states  0212,  0012  and  0210  are  larger  than  that  from  0010  is  indicative  of 

the  importance  of  rotational  transitions  in  vibrational  excitation.  Molecules 

in  the  states  where  and/or  j2^0  can  apparently  utilize  different  and 

9b 

stronger  portions  of  the  interaction  potential. 


B-  21  - 


Now  consider  the  general  V-R-T  relaxation  process  between  the  ground 
and  first  excited  vibrational  states 

H2(lji)  + H2(Oj2)  - H2 (0 j ) + H2  (0j2* ) 


The  rate  of  change  in  the  number  density, 
can  be  written  as 


of  molecular  species 


Mu,) 


dp 


Ik  _ 

at 


I 


- k 


kk'j; 


lji0j2-»0j  i'Ojj 


Pljt  P0j2  + k0j I'Ojj’  ljtOj2  Pojj'  P0j2’  (17) 


Since  rotational  relaxation  is  generally  much  faster  than  vibrational 

Relaxation,  it  is  reasonable  to  assume  that  the  rotational  states  form  a 
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Boltzmann  distribution.  We  may  therefore  write  the  number  density  as 


where  , . 

(23+1)  «p(-(E  -E  )/«) 

Pnj  J 

X (2J--+1)  exp(-  (Eni.-Eno)/KT) 

j' 


(18) 


is  the  rotational  popt lation  in  vibrational  level  n.  Substituting  Eq.  (18) 
in  Eq.  (17)  and  summing  over  leads  to 
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th 

dt 


V 1!ul 

L,  dt 
ji 


I 


all  indices 


' ^Oj^Oji'Oj,’  P1  Plj1  p0  P0j2 


+ ^ji'Ojj'^lj^jj  P0  P0JX*  P0j2*  (19) 


Equation  (19)  naturally  leads  to  a definition  of  the  effective  rate  constant  as 


r r 

,coi 00  = L Plji  P0J2  kljlOj2-»Oj1'Oj2( 


(20) 


The  effective  vibrational  equations  governing  the  process  are  then 


<*P, 


dt  ~ k01^00P0Pl  " k00-*0l  P0 


dPr  dp. 


dt 


where  the  second  equation  is  readily  established. 
This  set  of  equations  can  be  easily  solved  to  yield 


P0(t) 


k01->00  po^o) 


{[p0<°>  {k01--00  + k00-^0ll+  [ 


— P “ Pa(°)  {^01->00  * koo-oi}] 


oo-»oU  Loi->oor  rovv"  l 01  -> 00 


x e*p("  koi-»oopt)t 
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where  p = pQ+p^  *s  tke  tota^  density. 

It  can  readily  be  shown  that  as  the  system  nears  equilibrium 
(p0(0>  ~ Pq(®)  an<*  P^(0)  **•  P^C®))*  the  above  equations  can  be  linearized 

to  produce  the  following  simple  behavior 


pn(t)  - Pn(»)  = [pn(»)  - Pn(0)j  exp  (-  pQ(*>)  \t) 

for  n = 0,1.  The  relaxation  rate  in  the  exponent  under  these  conditions  is 
the  usual  sum  of  up  and  down  rates  in  a two-level  system 


1 ~ k0l->00  + koo->oi 

This  relaxation  rate  is  shovm  in  Fig.  7 labeled  vibrotor  (with  the  VB 

potential).  It  is  interesting  to  note  that  while  the  temperature  dependences 
tot 

of  k , , in  the  inset  of  Fig.  7 have  mostly  negative  curvatures,  the 

nlJln2J2 

averaging  procedure  in  Eq.  (19)  yields  an  overall  rate  with  the  normal 
positive  curvature.  The  vibrotor  calculation  clearly  lies  closer  to  the 
experiment  than  the  breathing  sphere  curves,  but  it  still  Falls  short  of 
precise  agreement. 
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VI.  CONCLUSION 

There  seems  to  oe  little  doubt  that  the  essence  of  the  relaxation 
mechanism  is  contained  in  the  full  vibrotor-vibrotor  treatment  with  the  VB 
potential.  Further  refinements  are  possible  with  more  sophisticated 
potentials  and  perhaps  improved  dynamical  methods.  An  important  conclusion 
from  the  above  analysis  is  that  rotations  do  play  a role  in  relaxation  (at 
least  in  the  low-temperature  region  examined  in  detail)  and  the  artificial 
exclusion  of  j ^ 0 states  in  the  breathing  sphere  approximation  under- 
estimates the  efficacy  of  vibrationally  active  collisions.  This  is  in  accord 

with  previous  work  on  He-H2  and  the  recent  parallel  semiclassical  treatment 
2c 

of  H2-H2.  Clearly .rotational  transitions  would  be  expected  to  continue  to 

enter  into  vibrational  inelasticity  at  higher  temperatures.  However,  this 

does  not  imply  that  the  breathing  sphere  rates  will  have  no  region  of 

5 6b  c 

applicability.  In  experiments  such  as  those  above  ’ * ' the  rotationally 

summed  and  averaged  rates  of  Eq.(l9)are  the  ones  of  relevance.  The  proper 
comparison  is  therefore  between  the  breathing  sphere  rates  and  those  of 

23 

Eq.(l9).  Model  calculations  have  shown  that  a region  of  applicability  exists 

for  this  approximation,  depending  on  the  nature  of  the  interactions.  Indeed 

22  25 

in  the  case  of  He-CO  * breathing  sphere  calculations  were  shown  to  be 

adequate.  Further  work  is  still  needed  on  this  important  problem. 

A few  final  comments  are  in  order.  These  calculations  are  moderately 

expensive  in  terms  of  computer  time  required.  It  should  be  recalled  that 

3b 

these  were  minimum  size  effective  Hamiltonian  calculations.  It  is  therefore 
easy  to  grasp  the  magnitude  of  difficulty  involved  in  studies  of  more  compJax 
(massive)  molecules.  Several  ways  of  circumventing  this  problem  can  be. 
suggested, ^ bu;  it  is  beyond  the  scope  of  this  paper  to  delve  into  this 
prob lera. 
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FIGURE  CAPTIONS 

1.  Body-fixed  axes  for  the  system  of  two  molecules.  R connects  the  centers 
of  mass  of  the  molecules  and  the  six  coordinates  are  R,  rl5  r2,  9j.>  0 2 
and  the  dihedral  angle  cp,  -y2. 

2.  (a)  Energy  spacing  for  the  breathing  sphere  levels  studied. 

- (b)  Energy  spacing  for  the  vibrotor  levels  considered. 

3.  (a)  Breathing  sphere  deexcitation  cross  sections  a(njn2  -» n1,n2' ) . Inset 

numbers  denote  the  quanta  lost  to  translation,  An  * (n^n-,)  - + n2* ). 

(b)  Selected  cross  sections  between  breathing  sphere  levels,  a(n>n2-» n^n^ ) . 

4.  (a)  Selected  rotor  cross  sections  in  both  ground  and  excited  vibrational 

levels. 

(b)  Profile  of  rotor  cross  sections  in  the  ground  vibrational  state  at  the 
total  energy  1.06  eV.  Lines  connect  the  same  initial  state,  0ji,0j2. 

5.  Profile  of  various  vibrotor  cross  sections  at  the  same  total  energy  1.2  eV. 
Lines  connect  the  same  initial  state  Oj1lj2. 

6.  Deexcitation  rates  for  H2  breathing  spheres.  Bands  encompass  rates  for 
transitions  involving  the  same  number  of  quanta  lost  to  translation, 

An  = (n,  + n2)  - (nt'  + n2’ ). 

7.  Rates  of  deexcitation  in  para-H2.  Results  of  calculation  with  the  London 
and  VB  potentials  are  compared  with  experiments . Vibrotor  results  are  only 
presented  for  the  VB  potential.  The  inset  shows  low  temperature 

summed  rates  defined  in  Eq.  (16). 
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FIGURE  3 (b) 
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ABSTRACT 


lb  this  wo rk  stochastic  theory  is  applied  to  the  treatment  of  atom- 
vibrotor  collisions.  This  is  an  extension  of  a previous  paper  which  described 
molecular  collisions  by  a Pauli  master  equation  or  a Fokker-Planck  equation. 

In  this  framework  an  energy  conserving  classical  path  model  is  explored,  and 
methods  for  solving  the  equations  r americally  are  discussed . The  coefficients 
of  the  Fokker-Planck  equation  are  shown  to  be  expressible  as  simple  functions 
of  the  interaction  potential.  Estimates  of  the  computational  labor  are  also 
discussed.  Finally  as  a follow-up  on  the  initicl  work,  numerical  solutions 

i 

of  the  master  equation  for  the  collinear  vibrational  excitation  problem  of 
Secrest  and  Johnson  are  presented  in  an  Appendix. 
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I.  INTRODUCTION 

Many  Important  and  Interesting  gas  phase  phenomena  involve  collisions 
of  relatively  large  and  complex  particles.  A previous  paper*  (hereafter 
referred  to  as  I)  investigated  treating  molecular  collisions  as  a process  of 
probability  diffusion  between  quantum  states.  This  avenue  of  approach  seems 
attractive  for  handling  complex  collision  systems  that  are  too  difficult  to 
treat  by  standard  methods.  These  ideas  are  more  fully  developed  in  the 
present  paper. 

The  treatment  In  I presented  equations  that  are  applicable,  in  principle, 
to  arbitrary  systems.  However,  the  theory  was  only  developed  in  detail  for 
or.e -dimensional  problems  (i.e . , one  internal  degree  of  freedom).  The  present 
paper  concentrates  on  an  atom-vibrotor  collision  system  as  the  simplest 
njeaningful  case  with  more  than  or.<*  degree  of  freedom.  The  resulting  stochastic 
equations  then  have  two  spatial  variables,  one  for  vibration  and  one  for  rota- 
tion. All  the  relevant  concepts  can  then  be  generalized  in  a straightforward 
fashion  to  arbitrary  inelastic  scattering  problems. 

The  stochastic  theory  requires  the  solution  of  either  a master  equation 
(ME)  or  a simpler  Fokker-pianck  equation  (FPE).  One  of  the  most  desirable 
aspects  of  this  theory  is  that  all  the  physics  of  the  collision  system  can  be 
condensed  into  a small  number  of  coefficient  functions  in  a FPE.  Qualitative 
predictions  about  the  behavior  of  the  inelastic  processes  can  then  be  imme- 
diately made  by  considering  the  magnitude  and  form  of  these  functions.  It 
will  be  shown  In  this  work  how  the  FPE  coefficients  can  be  generated  in  a 
straightforward  fashion  from  a given  intermolecular  potential.  This  approach 
may  ultimately  prove  to  be  useful  for  relating  observed  cross  sections  to  the 
properties  of  an  individual  collision  system  (i.e.,  the  Hamiltonian). 


In  most  cases  of  practical  interest,  the  ME  or  FPE  will  not  be  analyti- 
cally solvable.  This  paper  therefore  explores  methods  of  solving  these 
equations  numerically.  The  translational  degree  of  freedom  (represented  by 
the  vector  R)  is  treated  in  this  context  as  following  an  energy  conserving 
classical  trajectory.  For  most  problems  of  chemical  interest  the  particles 
are  massive  enough  so  that  the  classical  path  approximation  is  expected  to  be 
very  good.  The  combination  of  a classical  translational  path  with  internal 
degrees  of  freedom  obeying  quantum  stochastic  equations  of  motion  provides  an 
attractive  conceptual  model  for  molecular  collisions. 

Section  II  deals  with  the  detailed  consideration  of  the  stochastic 
theory  as  applied  to  atom-vibrotor  scattering.  It  is  shown  that  the  ME  can 
be  approximated  by  a FPE  which  describes  the  flow  of  probability  between 
energy  levels  rather  than  quantum  numbers.  This  has  important  advantages, 

i 

especially  for  case«  such  as  asymmetric  rotors  where  the  quantum  numbers  are 
poorly  defined.  Expressions  for  the  FPE  coefficients  are  developed  in  terms 
of  matrix  elements  of  simple  functions  of  the  interaction  potential.  Applica- 
tions of  various  approximate  methods  to  the  stochastic  equations  are  considered 
in  Section  III, 

' Numerical  methods  of  solution  are  treated  in  Section  IV.  Given  an 
interaction  potential,  an  energy  conserving  classical  path  model  for  the  ME 
requires  that  an  exponential  matrix  be  calculated  at  each  step  of  the  R 
integration.  If  n is  the  total  number  of  quantum  states,  this  procedure 
will  typically  require  matrix  multiplications  whose  computational  difficulty 
varies  as  n3.  The  FPE,  on  the  other  hand,  can  be  solved  in  principle  by 
methods  whose  difficulty  is  more  dependent  on  the  number  of  degrees  of  freedom, 
N,  than  the  number  of  states.  The  relative  computational  efforts  can  be  of 
considerable  practical  importance. 
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Finally,  calculations  were  presented  in  1 for  a collinear  vibrational 
excitation  problem.  It  is  shown  in  this  paper  that  some  of  the  approximations 
involved  in  those  calculations  can  be  eliminated  or  improved.  In  particular 
the  classical  path  formulation  of  Section  II  is  illustrated.  Results  for  the 
model  He-H2  collinear  system  of  Secrest  and  Johnson  are  presented  in  an 
Appendix.  In  addition  it  is  shown  in  Section  II  that  the  stochastic  equations 
remain  very  easy  to  solve,  even  for  noncollinear  collisions,  whenever  the 
intermolecular  potential  has  the  form 

V * A(R)  + B(R)C(r) 

where  r is  shorthand  for  all  the  internal  degrees  of  freedom.  The  present 
paper  further  develops  the  ideas  presented  in  I,  but  much  additional  work  is 
still  to  be  performed. 
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II.  STOCHASTIC  THEORY  OF  ATOM-VIBROTOR  COLLISIONS 

It  was  shown  in  I that  the  ME  for  a collision  system  with  L internal 
degrees  of  freedom  is  given  by 


P »X2 * * ’ * ,XN’ • * 


where  x^  Is  a (discrete)  index  for  the  i-th  degree  of  freedom.  If  the 
microscopic  transition  rates  A are  strongly  peaked  about  x^  * x^ , then 

P ^x^,x^, . . .,x^,t^  can  be  expanded  in  an  N dimensional  Taylor  series  about 

the  diagonal  values.  The  FPE  is  obtained  by  retaining  in  this  expansion  only 
terms  through  second  order.  It  will  therefore  contain  N terms  of  the  form 


N 


terms  of  the  form 


and  (N  - 1)  terms  of  the  form 


5*P 

3XjdXj  ’ 


Since  the  mixed  partial  derivative  cross  terms  only  occur  pairwise,  all  the 
essential  features  are  contained  in  a two-dimensional  example.  In  the 
remainder  of  this  paper  we  shall  use  an  atom-vibrotor  as  a prototype  two- 
dimensional  problem. 


A.  The  ME  for  an  Atom-Vibrotor  System 
The  relevant  ME  is 


Y.  A>J(n> j;n',j‘;t)  pj(n’,j',t) 

ri,  V 


atp 


(2.2) 


i 


$ 

* C - 6 - 

where  n is  the  vibrational  quantum  number,  j is  the  rotational  quantum 
number,  and  J Is  the  total  angular  momentum.  The  microscopic  transition 
rates  in  Eq.(2.2)  are  given  by  (see  I) 

i 

i 

r AJ(n,  j;n,.,j';t)  « J <n,  j j exp(-  iv  t/  ft)|  n',  j ')  | - 6 nn » 6 j j • (2*3) 

where  V is  the  intermolecular  potential  and  T is  the  (generally  time- 

dependent)  increment  betveen  time  steps.  Normally  the  orbital  angular  momentum 

i t should  be  included  in  Eq.(2.3),  but  it  is  assumed  to  be  eliminated  by  using 

2 

an  effective  Hamiltonian.  Note  that  V is  a function  of  the  translational 

— ^ ^ 

coordinate  R so  that  Eqs.(2,2)  and  (2.3)  assume  that  R is  a known  function 

of  time. 

It  was  shown  in  I that  the  above  ME  involves  approximations  whose 
validity  depends  on  the  size  of  t.  Most  important  of  these  is  the  use  of 
the  repeated  randomness  assumption  in  the  strong  interaction  region  of  the 
collision.  This  requires  that  the  probability  amplitudes  accumulate  enough 
phase  in  the  time  interval  t to  t + T for  the  random  phase  approximation 
to  be  made  at  each  step.  If  t is  too  small,  the  phase  accumulation  will 
not- be  large  enough  and  the  ME  will  not  be  a good  approximation.  Thic 
restriction  does  not  apply  to  the  weak  interaction  region  where  only  a small 
fraction  of  the  inelasticity  is  presumed  to  occur. 

There  are  also  constraints  on  the  largest  permissible  value  of  t.  The 
time  interval  t should,  clearly,  be  short  compared  to  the  duration  of  the 
collision.  It  may  also  be  noted  that  the  time  derivative  in  Eq.(2.2)  was 
obtained  as  a small  t approximation  to  the  finite  difference  expression 


T 


(2.4) 
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The  above  arguments  indicate  that  t aas  an  optimum  value  determined 
by  the  characteristics  of  the  system  at  each  point  in  time.  It  was  suggested 
in  I that  the  average  value  of  the  time  required  for  a jump  between  adjacent 
quantum  states  was  a good  estimate  for  t.  This  should  also  be  reasonably  close 
to  the  shortest  time  interval  for  which  the  difference  equations  are  numeri- 
cally stable.  If  Ho  is  the  Hamiltonian  for  the  internal  modes  (n  and  j) 
and  are  its  eigenvalues,  the  prescription  for  T(t)  becomes 

m>., 


il  (v'JlVrlf 

rf.j’ 

In  Eq.(2.5),  Ae^  should  be  interpreted  as  the  mean  of  the  energy  gaps 
between  nj  and  the  adjacent  quantum  states.  The  best  form  for  this  is  not 
entirely  clear  at  present  although  there  is  no  difficulty  for  a pure  vibra- 
tional problem  or  a pure  rotational  problem.  In  these  cases  the  geometric 
mean  of  the  two  nearest  energy  differences  is  indicated.  For  a harmonic 
oscillator  Asn  just  reduces  to  the  constant  energy  spacing  h u>.  The 
stability  criterion  of  Eq.(A.ll)  may  provide  a suitable  alternative  prescrip- 
tion for  t in  the  atom-vibrotor  case.  Practical  calculations  would  also  be 
very  helpful  in  this  regard. 

The  HE  is  still  not  complete  until  R(t)  is  known.  One  way  of 
specifying  this  is  to  assume  that  the  translational  degree  of  freedom  follows 
a classical  trajectory.  If  the  coupling  between  translation  and  internal 


T(t)  - £ P(n,J,t) 

n»  J 


- ft  £ P(n,j,t) 

n»  J 
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mode 8 is  then  ignored,  simple  models  can  be  developed  for  the  inelastic 
scattering.*-  A simplified  form  of  Eq.(2.5)  such  as 


T(t) 


ftA  c 


n0  jo 


[ Z (€n0jo 

n\j' 


(2.6) 


(n0,Jo  is  the  initial  state)  would  be  appropriate  in  this  case.  It  should  be 
noted  that  when  t is  an  appreciable  fraction  of  the  collision  time,  the 
time  derivative  in  Eq.(2.2)  should  properly  be  written  as  a finite  difference 
(see  Eq.(2.4)).  In  this  case  calculations  which  treat  (2.2)  as  a differential 
equation  in  time  may  significantly  overestimate  the  interaction  as  a result. 
However,  both  the  translational  decoupling  and  small  t approximations  are 
quite  good  in  the  limit  of  a large  number  of  strongly  coupled  internal 
states . 

Given  the  prescription  in  Eq.(2.6),  the  ME  is  very  easy  to  solve  if  V r 
is  time  independent.  It  is  easy  to  see  that  this  will  be  true  whenever  the 
potential  has  the  general  form 


V - A (t)  + B$)  C(r ) 


(2.7) 


where  r indicates  the  relative  separation  vector  of  the  diatom.  Inserting 
Eq. (2.7)  into  Eq.(2.6)  gives 


■ f)  A* 


noJc 


[ l (' 

n'.j' 


n0  jo 


' Vj')1  I' 


,2  "1-1/2  _k 

,|  J /B(R)  (2.8) 


n0 jo»n* j 


T (t ) 


Thus  the  rates  in  Eq.(2.3)  become 


T (t ) AJ(n,  j;n\ j';t) 


<n,j|exp 


Infill 


I ('no jo"  Vj")  fCn0j0)n«'j"l 
nH>  j" 


" 6nn'  6jjT 


All  of  the  time  dependence  in  the  rate  matrix  A is  now  seen  to  be  contained 

in  the  multiplicative  factor  l/r(t)  so  that  the  solution  of  Eq.(2.2)  is  the 

exponential  of  a time  independent  matrix  multiplied  by  a function  of  time. 

These  simple  models  have  the  disadvantage  that  the  classical  trajectory, 

$(t),  is  not  coupled  to  the  internal  degrees  of  freedom.  Since  this  will  not 

conserve  the  total  energy,  the  results  can  be  unsatisfactory  when  the  energy 

transfer  in  the  collision  is  an  appreciable  fraction  of  the  translational 

energy.  One  way  of  remedying  this  problem  is  to  force  conservation  of  energy 

3 

by  defining  the  instantaneous  radial  velocity  as 

• _ . 1/2 

— - ± {2/n  [e  - <i*)/2^R2  - <Ho>  - <V>]  } (2.10) 

Practical  application  of  this  approach  will,  of  course,  require  that  R(t) 

be  computed  numerically  from  Eq.(2.l0)  while  the  internal  modes  are  propagated 

dR 

by  Eq.(2.2).  The  sign  of  ^ is  not  determined  by  Eq.(2.l0)  so  that  it  may 
be  advantageous  to  integrate  the  momentum  conjugate  to  R as  well  (see 
Section  IV).  Expectation  values  of  Hq  and  V are  conveniently  given  by 


i 


10 


<Ho>  - £ «nj  P(n,j,t) 


(2.11«) 


<V>  - I £ p(n’:j',t)  <n,j|v(R)|n',J*>  (2.11b) 


n,j  rf.J1 


although  other  choices  are  possible.  The  orbital  angular  momentum  i may  be 

o 

treated  as  a constant  within  an  effective  Hamiltonian  formulation,  or  (l2) 
can  be  expressed  in  a proper  coupled  angular  momentum  representation. 


B.  The  FPE  for  an  At om-Vlb rotor  System 


For  this  special  case  the  FPE  is  obtained  by  a Taylor  expansion  of 
P~  (n ' , j ' , t ) in  Eq.(2.2).  However,  this  expansion  could  either  be  in  terms 
of  the  quantum  numbers  n,j  or  the  corresponding  energy  states  en,e^. 
can  be  shown  that  these  two  approaches  are  equivalent  in  that  their  finite 
difference  approximations  both  match  the  ME  through  second  order  (see 
Section  IV ).  However,  it  will  be  shown  below  that  if  the  expansion  is 
performed  in  terms  of  the  energy  levels,  the  FPE  coefficients  are  expressible 
as  matrix  elements  of  simple  functions  of  the  interaction  potential.  This  is 
in  contrast  to  the  development  in  I,  which  utilized  expansions  in  the  quantum 
numbers.  The  two  methods  are  identical  whenever  the  energy  levels  are  evenly 
spaced  as  for  a harmonic  oscillator. 
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The  FPE  in  energy  variables  is  now 


‘ B‘a^pJ(vVe)  + B*  ri^pJ(wt) 

+ '‘^(vy)  + CijfpJ(vvt)  + "ra-'lvj'*)  <2-12> 

J j n j 


where  the  coefficients  are  given  by 


■ kf  1 (*n*  * *n) 


(2.13a) 


kl  Z (C j * “ *J  ) AJ(n»JJn’»^,5t) 


(2.13b) 


I (en«~‘n)  AJ(n,j;nS 


n'.j’ 


(2.13c) 


The  sums  in  Eqs. (2,13a, b,c)  can  be  performed  analytically.  First  define  m 
to  be  the  reduced  mass  of  the  diatom  and  u ■ (r  - r ) to  be  the  vibrational 
displacement.  The  internal  Hamiltonian  Ho  can  be  decomposed  as 

H.  - H + H 
o v r 


Hv  “ pu/2m  + Vo(u) 


(2.14a) 


Hr  - j 2 /2  rar2 


(2.14b) 
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and 


H>.J>  - «>.!> 

(2.15a) 

Hr|n,J>  - 

(2.15b) 

These  last  two  equations  assume  that  rotation  and  vibration  are  separable, 
but  this  approximation  can  be  easily  removed.  The  various  s's  in  Eqs. (2.13a, b,c) 
can  now  be  replaced  by  and  Inserting  Eq.(2.3)  into  the  expression 

for  Bi,  Eq. (2.13a),  yields 


r'&i  » £ <n,j|exp(iVT/fi)ln,,J,><n,,j,|(Hv-en)  exp(-iVT/ft)|n,j>  (2.16) 

n',j' 

Noting  that  ^ |n', j ')(n*, j ' j is  just  the  unit  operator  leads  to  the  result 

n'.j' 


tBi  - <n,jiexp(iVT/ft)(Hv-en)  exp(-iVT/^i;!n,j)  (2.17) 

It  is  convenient  to  define  new  states  j^)  ana  jy  .)  by 

|fy>  « exp(iVT/ ft)(Hv  - <n)  exp(  - IVt/ fi)  | n,  j>  (2.18a) 

|tr^  - exp (iV T / fi)(Hr  " e j ) exp ( - IVt  / ft)j n,  j> 


(2.18b) 
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This  allows  Eqs. (2.13a, b,c)  to  be  compactly  written  as 


TBX 

- <n*j|vv> 

(2.19a) 

2tB2 

= <»J»v> 

(2.19b) 

TC: 

- <n,j|4rr> 

(2.19c) 

2t  C2 

“ <*r|*r> 

(2,19d) 

T D 

- <*rl*v> 

(2.19e) 

Furthermore  Eq.  (2.18a)  can  be  rewritten  as 

i 

Uv>  * exp(iVx/ft)  [(Hv-en),  exp(-iVr/fc)]  |n»j> 


where  [•,*]  is  the  commutator  bracket.  Therefore,  it  is  easy  to  see  that 


IV 


♦ T*(S?V  - 2t(^ 


\SmJ 


Su  / 


>u}  fn’j> 


(2.20a) 


and  in  a similar  fashion 


{t2  (1  - x2)  ( 


avy 

dx  / 


av  I 
dx  J 


av 

dx 


G 


2 t av  . , 

(i  - x2 ) a<p  •'z 


t2  / av\2 
(1  - x2 ) \ dcp  / 


i T ft 

(T^T 


(2.20b) 


+ 


where  x = cos  8 and  the  operator  G is  defined  by 


G|j>  - Jiiiil  {ilc» 

/2j  + 1 V2j  - 1 


j + l>  \ 


2j  + 3 


The  three  ep  dependent  terms  in  Eq. (2.20b)  can  be  removed  by  transforming 

to  the  body-fixed  (BF)  reference  frame^  or  by  using  at*  ef fee;  Ive  Hamiltonian 
2 

approximation. 

The  new  states  jty  ),  |^r)  have  an  interesting  physical  interpretation. 
From  the  Heisenberg  , equations  of  motion 


r ft  LSx 


Sx  (l-x2)  Sep2  (1-x2)  Sep  Jz 


The- states  ♦ ) can,  therefore,  be  expressed  in  the  form 


IV-  K + £(!*)’] 


2mr2|*t>  - [ri(  + (§'/  + T!(l-x*)  (g)!]|n,J>  (2. 


which  are  of  the  form  of  t times  a rate  of  energy  transfer  plus  t2  times 
a dispersive  term. 


III.  APPROXIMATION  METHODS  AND  THE  STOCHASTIC  THEORY 


Various  standard  approximation  schemes  for  inelastic  scattering  can  be 

incorporated  in  the  stochastic  theory.  In  this  regard  it  has  been  shown  above 

that  the  FPE  coefficients  can  be  obtained  directly  from  the  potential.  The 

properties  of  these  coefficients,  without  further  computation,  may  then  be 

used  as  a simple  guide  to  the  behavior  of  the  collision  system  or  may  provide 

physical  insight  into  the  approximation  methods.  For  example,  in  this  section 

we  shall  consider  the  combination  of  breathing  sphere  and  effective  Hamiltonian 

methods  with  the  stochastic  theory. 

Effective  Hamiltonian  methods  reduce  the  dimensionality  of  a close 

coupled  calculation  by  preaveraging  or  eliminating  angular  momentum  projections 

2 

before  doing  the  collision  dynamics.  The  centrifugal  decoupling  (CD)  approxi- 
mation^ is  easily  incorporated  into  a formulation  of  the  ME  or  FPE.  Use  of 

g 

an  effective  potential  with  the  FPE  is  inconvenient  unless  a operator 

can  be  defined,  thus  permitting  the  use  of  Eqs . (2.20a, b) . However,  it  has 

recently  been  shown  that  modified  effective  potential  operators  can  be  defined 

9 

which  retain  the  same  dimensionality  reduction  as  the  original  formulation. 

, In  the  BF  reference  frame  the  potential  has  no  dependence  on  azimuthal 
angles,  Eq. (2.20b)  therefore  simplifies  to 


2”r!|*r>  ” + tT*[J7<1"'S>s]  + 2lTft  tfGl  |n>1>  (3 


where  x is  now  the  cosine  of  the  angle  between  it  and  r.  The  quantum 
number  w (uhe  projection  of  the  rotational  ingular  momentum  on  the  BF  z 
axis)  should  also  be  included  in  the  ME,  Eq.(2.2),  in  this  case.  Then  (X2) 


ir.  Eq.(2.lO)  is  given  by 


<i*>  - J2  + <j2)  - 2 <iu2>  (3.2 

if  the  CD  approximation  is  used.  In  an  effective  potential  formulation,  on 
the  other  hand,  l2  is  treated  as  a constant  throughout  the  collision. 

The  breathing  sphere  approximation  assumes  that  there  is  no  strong 
vibration-rotation  coupling.  Defining  P ( *n» fc)  as 

f(ve)  - 1 Kvv') 

j 

it  is  easily  seen  that  p(en»t)  will  satisfy  a breathing  sphere  FPE  if  the 
vibrational  and  rotational  parts  of  Eq.{2.12)  are  separable.  This  will  be 
true,  for  instance,  if  D«0  and  if  BitB2  have  no  rotational  dependence. 
Since  the  coefficients  can  be  obtained  relatively  easily,  it  should  be 
possible  to  determine  if  a breathing  sphere  FPE  is  a good  approximation  to 
Eq. (2,12), 
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IV.  NUMERICAL  METHODS  FOR  SOLVING  THE  STOCHASTIC  EQUATIONS 

There  will  obviously  be  only  a small  class  of  scattering  problems  for 

which  the  ME  and  FPE  are  analytically  solvable.  Numerical  methods  for  solving 

the  equations  with  appropriate  boundary  conditions  are,  therefore,  considered 

in  this  section.  Consistent  with  the  goal  of  describing  the  scattering  process 

by  a simple  equation,  we  shall  assume  that  an  effective  potential  will  be 
8 9 

used.  * The  translational  degree  of  freedom  will  also  be  described  by  an 
energy-conserving  classical  path  as  outlined  in  Section  II. A. 


A.  Numerical  Methods  for  the  ME 


Within  the  classical  path  approximation,  the  system  is  started  in  an 
initial  state  nQJ0  at  a value  of  R outside  the  range  of  the  potential. 
For  each  value  of  the  orbital  angular  momentum  l,  the  equations  for  the 
Internal  modes  and  R(t)  are  integrated  until  the  particles  are  separating 
and  R is  again  outside  the  range  of  the  potential.  Each  such  "trajectory” 
will  yield  the  complete  vector  of  transition  probabilities  into  all  the 
possible  final  states  n^ 


Total  cross  sections  for  translational  energy  E1  are  then  computed  by 


n j -*n  j 
oJo  1J1 


" 2|iEi(2j0  + l)  I (21+1)  Pnojo-»niji(E)  & l) 
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Internal  state  probabilities  are  propagated  from  time  t to  t + T by 

J 

the  finite  difference  form  of  the  ME 

i 

t 

i 

i 

P^(n,j,t+r)  - l j<n,  jjexpt- iV(R)  tM]  } r\* , j *>  | P^(n',j',t)  (4.2) 

n’.j’ 

where  T is  given  by  Eq.(2.5).  The  translational  degree  of  freedom  is 
handled  by  numerically  integrating  the  equations 

R * sign(pR)  x{2[e  - fl2jfc(l  -i- 1)/2 p.R*  - (h0>  - (v) ] / p.j * ' * (4.3) 

pR  = (J.R  * h2i(l+  1)/ p,R3  - ^ (v)  (4.4) 

with  stqp  size  T.  Since  pR  is  only  needed  to  determine  the  sign  of  R, 

it  does  not  matter  that  Eq.(4.4)  will  not  conserve  energy.  At  the  end  of 

each  time  step  pR  can  be  set  equal  to  p.  times  R computed  from  Eq.(4.3). 

An  application  of  these  ideas  to  a collinear  vibrator  is  presented  in  the 
Appendix. 

When  the  dimensionality  n of  the  problem  is  reasonably  large,  the 
major  computational  expense  of  this  procedure  is  the  evaluation  of  the 
microscopic  transition  rates  on  the  R.H.S.  of  Eq.(4.2).  This  would  normally 
involve  the  numerical  exponentiation  of  a matrix  whose  elements  are 

- iT<n,j|v(R)|n,,j')/ft 


When  the  potential  has  the  form  of  Eq.(2.7),  the  time-consuming  matrix 
diagonalization  needs  to  be  performed  only  once.  For  a more  general  potential 
Eq. (4.2)  will  require  matrix  manipulations  whose  difficulty  varies  like  n3 . 
Computational  expense  can  therefore  be  a serious  problem  for  large  n,  but 
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B.  Numerical  Methods  for  the  FPE 

Solution  of  the  FPE  with  a classical  path  assumption  should  be  very 
similar  to  the  procedure  outlined  above  for  the  ME.  Internal  state 
probabilities  are  propagated,  however,  by  a finite  difference  form  of  the  FPE 
rather  than  by  Eq.(4.2).  The  FPE  may  be  viewed  as  an  approximation  to  the  ME 
in  the  limit  that  the  discrete  quantum  states  are  closely  spaced  enough  to 
resemble  a continuum.  The  possibly  oscillatory  transition  rate  matrix  A is 
also  approximated  by  a smooth,  continuous  distribution,  which  is  determined 
by  the  first  and  second  moments  of  A about  the  diagonal.  Consequently,  a 
numerical  solution  of  the  FPE  does  not  require  the  consideration  of  every 
state,  and  the  computational  difficulty  depends  more  on  the  number  of  degrees 
of  freedom  than  the  number  of  states.  This  can  be  a great  advantage  for 
complex  collision  systems  where  the  number  of  open  channels  Is  very  large. 

The  FPE  is  to  be  solved  on  a region  with  a stepped  boundary  since 
conservation  of  the  total  energy  links  the  maximum  allowed  vibrational  and 
rotational  levels.  This  area  must  then  be  filled  with  a mesh  or  grid  of 
discrete  points  at  which  the  partial  derivatives  of  the  FPE  ere  replaced  by 
finite  differences.  It  is  also  desirable  to  express  the  boundary  conditions 
in  a numerical  fashion  that  is  independent  of  the  particular  functional  form 
of  the  coefficients.  Since  the  FPE  is  derived  from  the  ME,  formulation  of 
appropriate  difference  expressions  and  boundary  conditions  may  be  guided  by  a 
comparison  with  the  ME.  It  is,  therefore,  suggestive  that  the  mesh  be  no 
more  closely  spaced  than  the  discrete  quantum  states.  A typical  (but 
simplified)  situation  is  shown  in  Figure  1 where  each  grid  point  is  assumed 
to  lie  on  one  of  the  quantum  states.  This  figure  is  analyzed  below. 
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Define  x and  y to  be  the  continuous  variables  corresponding  to  the 
vibrational  energy  and  rotational  energy,  respectively.  Now  consider  the 
state  11  of  Figure  1 and  the  eignt  states  surrounding  it  to  be  typical  of  a 
point  away  from  the  boundaries.  Since  the  FPE  was  obtained  by  an  expansion 
through  second  order  of  the  ME,  it  is  natural  that  the  continuous  function 
p(x,y)  be  obtained  by  quadratic  interpolation.  The  expressions  for  the 
partial  derivatives  are,  therefore,  obtained  from  the  biquadratic  Lagrange 
interpolating  polynomial*0 


(x-x  )(x-x  ) 

»y)  “ * 


xiox2o 


p00(y-yi)(y -yj)  _ p01(y-y0)(y-y2)  | P02(yy0)(y-yl) 
yloyio  ~ yioy2i  y2<>y2i 
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X1QX21 


P10(y-yx)(y  -y2)  ptl(y-y0)(y -y2)  p12(y  -y0)(y  -yx) 


yioy«o 


y,Dy 


lO'll 


y20y2i 


(x-xo)(x-xi) 
X2  0X2  1 


pa>(y  •y1)(y-y1)  p21(y-y0)(y-y2)  ?22(y  - y0)(y  ■ yx) 

yi0y2i  + 


yi0y 


l <K2  0 


y2ny 


(A. 5 


2oJ  21 


where  P^  « P^.y^)  and  x^  « xA  - x ^ . 

The  following  expressions  are  thus  obtained  by  differentiating  Eq.(4.5): 


3P 

X10X20X2l  55 


=**  V p 

X = Xl  X10r2l 


x2  P 
X21r01 


x20(xio  ' x71)  pn  (4-6a) 


y *y 
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X»X1 
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^ [_X2 1 P0 1 " X20Pll  + XlOP2lJ 


(4.6b ) 
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X * Xj  **  Xl0yi0P22  + X2iy21P00  ’ X2iyi0P02  ” X10y21P20 

ymyl 

+ X2iy20^yio"y21  )P01  + y21X2o(XlO  " X21  )P10 
Xl0y20  ^yiO  " y2  1 ^2  1 " ytOX20  ^X10  ' X21  ^12 
+ X20y20(Xl0  -X2l)(yiO  “y2l)PH  (4‘?) 

ui.1  g2p 

Approximations  for  ^ and  are  obtained  by  switching  x and  y in 

Eqs . (4.6a,b) . If  it  is  assumed  that  the  microscopic  transition  rates 
A(n, j; n* , j ' ; t)  are  zero  for  more  than  one  step  off  the  diagonal,  it  can  bs 
shown  that  the  finite  difference  form  of  the  FPE  is  then  identical  to  the 
appropriate  ME  in  Eq.(4,8) 


d2P 

:i0X20X2iyi0y20y21  dx&y 


^P(n,j,t)  - A(n,j;n+ 1, j,t)  [p(n+ l,  j,t)  - P(n,j,t)J 

+ A(n,  j;n  - 1,  j,t)  jj(n  - 1,  j,t)  - P(n,  j,t)] 

+ A(n,  j;n,  j+  l;t)  £p(n,  j+  l,t)  - P(n,  j,t)J 
+ A(n, j;n,J  - 1; t ) [p(n, j - l,t)  - P(n,j,t)]  (4.8) 
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In  this  connection  note  that  the  D coefficient  involves  microscopic  transi- 
tion rates  that  are  at  least  two  steps  off  the  diagonal  (one  in  n and  one  in  j). 

The  finite  difference  form  of  the  FPE  presented  here  therefore  matches 
the  ME  at  n-1,  n,  and  n+ 1 (and  similarly  for  j).  When  n is  equal  to 
zero,  modified  approximations  to  the  partial  derivatives  which  match  the  ME  at 
n,  n+  1,  and  n+  2 should  be  used.  Appropriate  revised  expressions  for 


dP  4 
- — end 
Ox 


the  left 


d2P 

bxdy 

edge 


are  obtained  by  evaluating  the  derivatives  of  Eq.(4.5)  along 
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e case  j»0  is  treated  by  analog,  _ ’d  at  the  corner  point  0 0 
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Suitable  boundary  conditions  also  need  to  be  formulated  for  the  upper 
(stepped)  edge.  These  must  satisfy  the  requirements  that  all  the  probabilities 
sum  to  unity  while  the  flux  ^ ~ or  ^ vanishes  everywhere  along  the 
upper  border.  This  last  condition  can  be  fulfilled  by  assigning  fictitious 
probabilities  to  the  layer  of  states  immediately  outside  the  region  of 
interest.  It  should  be  pointed  out,  however,  that  these  non-zero  closed 
channel  probabilities  are  only  a mathematical  construct  with  no  physical 
meaning.  Referring  to  Figure  1,  we  thus  set  P0,=P08,  Pl9«P18,  P29“Pi8’ 

P28  = pi 8’  etc*  a^ove  boundary  conditions  will  still  not  insure  conserva- 

tion of  probability;  this  can  be  accomplished  by  renormalizing  the  individual 
probabilities  at  each  step  in  time.  It  is  easy  to  see  that  the  procedure 
outlined  here  is  consistent  with  a statistical  solution,  i,.e., 

t 

p(n,j,t)  = 1/ (number  of  open  states) 

Leakage  (dissociation)  can  also  be  allowed  if  the  fluxes  through  part  of  the 
upper  boundary  are  different  from  zero  (see  paper  I). 

Reference  to  Figure  l also  reveals  that  the  zero  flux  condition  requires 

^12  an<*  Pjj  "P44  = P53  ~P62*  s^ou^  recalled  that  only  open  states 

are  included  in  these  calculations,  and  probabilities  along  the  upper  boundary 
may,  therefore,  be  unreliable.  This  is  a common  situation  with  any  finite 
basis  set  calculation.  In  the  present  case  no  difficulty  should  be 
encountered  since  the  stochastic  theory  is  developed  for  large  systems  where 
the  upper  border  is  likely  to  be  far  from  the  region  of  interest.  In  other  words, 
most  problems  will  be  concerned  with  behavior  in  the  lower  left  side  of 
Figure  1 and  the  presence  of  the  upper  boundary  will  frequently  not  be  sensed. 
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In  view  of  the  above,  it  is  easily  seen  that  propagating  in  time  the 
probability  of  being  in  any  specific  quantum  state  by  means  of  the  FPE  requires 
only  the  surrounding  block  of  nine  states.  Thus  the  numerical  effort  of 
computing  each  "trajectory"  for  this  two-dimensional  problem  goes  like  9n. 

More  generally  the  dependence  would  be  n X 3 (n  is  the  number  of  states 
and  N is  the  number  of  degrees  of  freedom).  Numerical  computation  of  the 
FPE,  like  classical  mechanical  methods,  depends  more  on  the  number  of  degrees 
of  freedom  than  the  number  of  quantum  states . This  should  be  a tremendous 
advantage  when  n is  very  large,  as  it  is  for  many  problems  of  interest.  In 
the  limit  of  many  states,  the  FPE  may  also  be  treated  as  an  exact  partial 
differential  equation  rather  than  as  a smoothed  approximation  to  the  ME.  The 
FPE  can  thus  be  solved  with  a coarser  mesh,  resulting  in  additional  savings  in 
computational  labor. 

Considerations  bearing  on  the  choice  of  the  time  steps  t were  discussed 
in  Section  II. A;  for  the  FPE  the  numerical  stability  of  the  difference  equations 
should  also  be  taken  into  account.  We  wish  to  apply  the  von  Neumann  stability 
criterion^  to  the  FPE.  Assume  P^xn»ym,t^  ^as  ^orm 

P(Xn,vVt)  = exP(in9)  exp(iracp)  exp(i\t) 

where  9,  cp , and  \ are  some  numbers.  The  difference  equations  are  then 
presumed  stable  if  exp(iXt)  is  a decreasing  exponential  for  all  real  9 
and  <p.  Applying  this  to  the  atom-vibrotor  FPE  yields  the  following 
inequality:  ^ 

{ [,X20  ^X2  l " X10  ^ 1 " ^ X20  TBj  J ^10^21 

+ [y20(y2i  -yio) rCi  - 2y2oTC2 J xu>x2t  + 2(*2i  -*io)(y2i  -yi0)TDj 


(4.11) 


where  the  indices  0,  1,  2 stand  for  n-1,  n,  n + 1,  and  similar  indexing  applies 
to  j.  The  maximum  value  of  t allowed  by  Eq.(k.ll)  is  thus  a simple  function 
of  the  FPE  coefficients . 


If  the  interaction  is  weak,  Eq.(4.11)  can  be  readily  compared  with  the 
previous  prescription.  Expanding  the  exponential  time  propagator  and  assuming 
a tridiagonal  potential  matrix  gives 


TA(n,j;n',j';t)  - j t21  j •!  / 


1 n - n'  | <1 
and  | j - j ' | < 1 


Inserting  this  into  Eq.(4.11)  yields 


otherwise 


T - lVnj,nH  jl  + lVnj,n->  jl  + lVnj,nj+J  + fVnj,nj-J  ] 


2 -1-1/2 


which  is  similar  to  the  physical  criterion  presented  in  Eq.(2.5). 
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V.  CONCLUSION 

In  this  paper  we  have  formulated  additional  practical  procedures  for 
dealing  with  the  stochastic  equations  presented  in  I.  Molecular  collisions 
frequently  involve  relatively  massive  particles  so  that  the  translational 
degree  of  freedom  may  be  reasonably  treated  as  following  a classical  trajectory. 
We  have  presented  a formalism  in  which  the  translational  degree  of  freedom  is 
described  by  classical  mechanics  while  the  internal  states  are  propagated  by 
a quantum  stochastic  equation  of  motion.  Considering  the  simplicity  of  the 
stochastic  equations  for  the  internal  modes,  this  would  seem  to  be  an 
attractive  conceptual  model  for  molecular  collisions.  Indeed,  useful  qualita- 
tive information  can  be  gained  from  a knowledge  of  the  FPE  coefficients,  even 
without  solving  the  equation,  since  these  coefficients  control  the  effective 
cpllisional  coupling. 

There  are  numerous  techniques  for  calculating  collisional  information, 
and  each  has  its  own  realm  of  applicability.  Previous  methods  have  generally 
proved  inadequate  for  dealing  with  the  scattering  of  relatively  complex 
particles.  For  many  systems  of  interest,  the  number  of  open  channels  n is 
too  large  for  a fully  quantal  close  coupling  or  effective  Hamiltonian  calcula- 
tion to  be  practical.  The  expense  of  performing  completely  classical  trajec- 
tory studies  depends  more  on  the  number  of  internal  degrees  of  freedom  N 

than  the  number  of  states.  They  have  the  serious  disadvantage,  however,  that 

13 

all  quantum  effects  are  lost.  Classical  S~matrix  methods  combine  quanftm 

interference  effects  with  the  ease  of  computing  classical  trajectories.  Even 

9 

with  the  use  of  effective  Hamiltonians,  there  may  still  be  too  many  degrees 
of  freedom  for  a classical  S-raatrix  calculation  to  be  feasible  for  many 
collision  problems.  It  is  hoped  that  the  FPE,  as  formulated  in  this  work, 
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will  permit  reasonable  and  practical  scattering  calculations  for  very  complex 
collision  systems. 

Finally,  the  FPE  is  a type  of  partial  differential  equation  which  has 
been  extensively  studied.  It  was  shown  in  I how  a knowledge  of  the  FPE 
coefficients  and  the  associated  boundary  conditions  could  lead  to  reasonable 
qualitative  predictions  about  the  behavior  of  an  inelastic  collision.  This 
paper  presents  a practical  nethod  for  proceeding  directly  from  the  interaction 
potential  tc  the  coefficients  of  the  FPE.  These  could  then  be  used  for  quali- 
tative interpretation  or  as  input  into  numerical  calculations.  It  is  therefore 
hoped  that  a reasonably  direct  connection  can  be  made  between  the  Hamiltonian 
of  a system  and  the  outcome  of  an  Inelastic  collision. 


ft 

1 
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APPENDIX:  CALCULATIONS  FOR  THE  COLLINEAR  VIBRATIONAL 

EXCITATION  OF  A MODEL  Re-H2  SYSTEM 


The  collinear  collision  of  an  atom  with  a diatomic  harmonic  oscillator 

has  been  used  frequently  as  a test  case  for  various  approximate  methods. 

Calculations  including  certain  simplifying  assumptions  were  presented  in  I. 

When  these  additional  approximations  are  removed,  the  solutions  of  the  ME  can 

14 

be  directly  compared  to  '’exact"  quantum  results,  and  this  is  done  below. 

The  collision  system  was  described  in  I;  as  before  the  interaction 
potential  is  taken  to  be  an  exponential  repulsion  between  the  atom  and  the 
near  end  of  the  diatom.  In  the  reduced  units  of  Ref.  14a  (except  that  the 
unit  of  energy  is  ftu)  not  h ta/2),  the  Hamiltonian  is 


E = H/ftu)  * P*/2m  + Hn/flu)  + V0  exp(-cvR)  exp(o-x) 


(4,1) 


where  HQ  is  the  vibrational  energy  of  the  oscillator,  x is  the  reduced 
displacement  of  the  oscillator  from  equilibrium,  R is  the  reduced  transla- 
tional coordinate,  and  P is  its  conjugate  momentum.  Since  the  Hamiltonian 

K 

14  b 

is  invariant  to  the  transformation 


R + 6 


VQ  may  be  arbitrarily  set  equal  to  E. 

The  finite  difference  form  of  the  ME  is 


Pn(t  + T)  = £ A(n,m,t) 


ra 
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where  t Is  in  units  of  u)'1  and  A'(n,m,t)  ^ TA(n,m,t)  + 6^  with 

2 

A^n,m,t)  = | (nj  exp[iTE  exp(-aR)  exp(o'x)]  |ra)  | (A. 3) 

The  matrix  elements  of  exp(orx)  are 

(n|exp(a  x)|m)  * (- 7)  (a  2"1/2 )m  n exp(a2/4)  L (-a2/2)  (A. 4) 

for  m > n (Ln(x)  is  the  associated  Laguerre  polynomial;.  Given  values  of 

T and  F.,  the  rates  in  Eq.(A,3)  are  then  conveniently  obtained  by 
numerically  exponentiating  the  matrix  whose  n,m-th  element  is 

iTE  exp(-O'R)  (n|  exp  (ax)  |m) 

The  ME  was  solved  by  the  procedure  outlined  in  Section  IV. A.  A fifth 
order,  varifble  step  size  Adams-Moulton  predictor-corrector  algorithm^  was 
used  to  compute  the  classical  trajectory  for  the  R degree  of  freedom.  At 
each  step  energy  was  forced  to  be  conserved  by  setting 

R(t)  - sign  (PR)  X {2^E  - < Hq / fi u>)  - E exp(-aR)  (exp(ax))  J /mj  ^ (A. 5) 

while  P,,  was  predicted  by  integrating 
K 

PR  = H-aEexp(-aR)  (exp(ax))  (A. 6) 
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from  t to  t + T (cf.  Eqs , (4.3),  (4.4) ) . The  expectation  values  in  these 
equations  are  given  by 

<H0/ftu>>  = ^ (n+l/2)Pn(t) 

n 

(exp (ax)  ) = l 7Pn(t)  exp(ofx) | m> 

n,m 


Since  P„  is  only  needed  for  determining  the  sign  of  R,  it  does  not  matter 
that  repeated  application  of  Eq.(A.6)  will  not  conserve  the  total  energy.  At 
the  beginning  and  end  of  each  step,  P /m  was  set  equal  to  R computed  from 
Eq. (A.5 ) . 

> Each  trajectory  is  started  in  some  initial  state  nQ  with  R set  equal 
to  some  large  value  outside  the  range  of  the  potential.  The  time  steps  t 
should  be  given  by  (cf.  Eq.(2.5)) 

T ( t ) ~ y Pn(t)  exp  ^-Q'R(t)J  ^(n-m)2  | (n  exp  (a  x)  |m)  | j j (A.  7) 

n m 

when  in  the  strong  interaction  region.  This  equation  is  not  applicable  for 

large  R where  the  interaction  is  weak  and  it  predicts  steps  much  too  large. 

As  a practical  matter  then  the  time  step  at  the  start  of  the  trajectory  was 

given  some  value  t , just  larger  than  the  estimated  minimum  time  from 

Eq.(A.7).  This  was  done  to  insure*  that  the  strong  interaction  region  was  not 

bypassed  in  the  integration.  When  Eq.(A.7)  predicted  a smaller  value  than 

T , this  smaller  value  was  used, 
o 
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Calculations  were  performed^  for  the  parameter  set  m=2/3,  O' = 0.314 
at  E = 8.0;  and  for  m=2/3,  cv  = 0.30  at  E = l0.0.  The  two  values  of  a 
were  necessary  because  of  differences  ii  the  calculations  of  Refs.  (14a)  and 
(14b).  Results  are  plotted  in  Figures  2 and  3 along  with  the  corresponding 
quantum  values.  It  may  be  seen  that  the  stochastic  curves  seem  to  represent 
the  quantum  mechanical  transition  probabilities  with  the  oscillations  averaged 
out.  This  result  is  not  surprising  in  view  of  the  removal  of  phase  interfer- 
ences in  the  stochastic  theory.  It  should  also  be  pointed  out  that  the  col- 
linear  harmonic  oscillator  model  apparently  exhibits  an  anomalously  large 

amount  of  oscillation.  Indeed  quantum  calculation  for  a Morse  oscillator 

Ikb. 

show  much  less  of  this  kind  of  structure. 

Calculations  were  also  carried  out  at  energies  below  E = 8.  However, 
these  results  were  in  generally  poor  agreement  with  the  quantum  values  for  E 
less  than  6.0.  This  is  not  unexpected  since  the  stochastic  theory  should  be 
best  where  there  are  a large  number  of  strongly  coupled  states . 

The  procedure  used  here  for  conserving  energy  no  longer  guarantees  the. 
satisfaction  of  microscopic  reversibility  (i.e.,  _^n  (E)  - pn  (E)j. 

However  f.:~  classically  allowed  transitions,  the  probabilities  for  excitation 
and  deexcitation  were  still  in  reasonable  agreement.  In  contrast,  the 
probability  of  deexcitation  is  much  smaller  (and  more  realistic  than  the 
excitation  probabilities)  when  the  transitions  were  classically  forbidden. 

The  origin  of  this  phenomenon  is  unclear. 
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Figure  Captions 


Figure  1. 

A set  of  mesh  points  for  the  solution  of  an  atora-vibrator  FPE.  The 
equation  is  to  be  solved  in  the  shaded  region  while  the  layer  of  closed 
states  outside  of  the  upper  edge  is  used  to  formulate  the  boundary  condi- 
tions. The  number  of  states  and  their  spacings  are  intended  for  illustra- 
tive purposes  only  and  should  not  be  construed  to  represent  any  real  system. 
The  continuous  variables  x and  y correspond  to  the  vibrational  energy  and 
rotational  energy,  respectively. 

Figure  2. 

Plot  of  transition  probabilities  from  initial  states  no  as  a function 
of  final  state  n (Pn  ) for  E = 8.  These  are  for  the  collinear  oscillator 
problem  with  a = 0.314  and  m = 2/3-  The  solid  lines  connect  the  quantum 
values  of  Clark  and  Dickinson*^*5  while  the  open  circles  are  the  stochastic 
results  of  this  work. 

Figure  3- 

The  same  as  Figure  2 for  E = 10  except  that  a = 0.30.  Quantum  results 

l4a . 


are  those  of  Secrest  and  Johnson 


